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SECTION I 
INTRODUCTION 


The concept of supporting a rapidly spinning body by means of electric 
fields to obtain an inertial reference is motivated by the desire to avoid 
a common problem among all conventional inertial devices - mechanical 
friction that occurs among the linkages that support the spinning body (or 
bodies). Fluid-bearing devices have been used as a partial answer to the 
problem and serve well in short-term applications in such severe environ- 
ments that the use of conventional sensors is impractical, if not impossible. 
Viscous forces associated with these devices generally render them un- 
suitable for long-term applications. 

By supporting a spherical body (the rotor) by electric fields in an evacuated 
housing, mechanical and viscous forces on the rotor are removed. Long- 
term stability of the momentum vector, or compensated momentum vector, 
becomes adequate for many long-term inertial reference and navigation 
applications. However, mere establishment of such a stable inertial ref- 
erence is not adequate; to gainfully utilize this reference, one must be able 
to measure exactly where that vector is, relative to a convenient set of 
coordinates. 

The ideal mechanism for measuring the momentum vector of the rotor 
(reading out the gyro) is one which does not degrade the stability and pre- 
dictability of the device. In a strict sense, this mechanism does not exist; 
therefore, the best alternate is one which minimizes the disturbance. An 
optical device is one such mechanism in that the minimum disturbance is 
the impulse of photons that strike the rotor; this disturbance is of the order 
of the ratio of the photon momenta to the rotor momentum, which is much 
smaller than the disturbance offered by the torque produced by electric 
supporting fields. 

To employ an optical system, the rotor must possess some means of 
codifying the optical signals received by the device; the information can 
then be processed to determine the direction of the momentum vector. A 
spinning rotor automatically "codes" the signal by a doppler shift of fre- 
quency of the incident signal by the moving surface. Measurement of this 
shift determines the location of the spin vector; if it coincides with the 
momentum vector, readout is accomplished. However, the technology of 
measuring the doppler shift to the required accuracy has not advanced to 
the point of feasibility for the ESVG. Therefore, another method of 
codifying the signal is needed - at least for the present time. 
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The present method of codifying the optical signal is to manufacture the 
rotor so that its surface is highly reflective in some areas and diffuse 
and/or absorptive in other areas. Electric torque and electric field 
breakdown considerations dictate that the diffuse and absorptive areas 
(where the surface finish is rough and jagged in a microscopic sense or 
where material is deposited that exhibits different electrical properties) 
be kept to a minimum. The diffuse areas define a pattern in such a 
manner that the code generated by it can be related uniquely to a specific 
colatitude of the surface with respect to the pattern axis. This pattern 
axis is aligned as closely as possible to the momentum axis first by 
fabricating a preferred axis (axis of maximum moment) into the rotor, 
then by aligning the pattern axis to be coincident with the preferred axis, 
and finally by damping out the initial rotor nutation in the startup process. 

Two approaches of signal codification have been employed on strapdown ESVG 
readout systems. One approach uses coded latitude lines; a specific sequence 
of pulses on a parallel of latitude is used as a coarse estimate of the region 
of the rotor on which the optical signal is focused. A series of alternating 
lines provides a fine measurement by measuring the amplitude of the 
signal. This type of pattern occupies a considerable area of rotor surface 
and fails to minimize the diffuse areas of the surface. 

Another approach utilizes fine pattern lines in which the pattern edge precisely 
defines a specific function of the colatitude. The colatitude is uniquely 
determined by the longitude angle that separates the points of two pattern 
lines viewed by the optical system and the functions which the pattern lines 
define. 

It is the latter type that is presently employed on a number of strapdown 
ESVG's. There is no evidence in recent test data from these units that 
this type of pattern line, as presently applied, degrades the ESVG drift 
performance. However, the same test data indicate that the accuracy of 
the readout system needs to be improved to be commensurate with the 
predictability of the ESVG momentum vector and to match accuracies of 
readouts on other inertial devices. The need for improvement, plus the 
need to streamline the information processing functions so that it may be 
implemented on a vehicle, is the motivation for this study program on the 
strapdown ESVG readout system. 

The program was based on the premise that the approach presently em- 
ployed (optical system, longitude measurement between pattern lines) has 
the capability of reading out the ESVG momentum vector direction to an 
accuracy of the one-second range. It was first necessary to identify the 
dominant sources of error, which if reduced would significantly improve 
system error. An error analysis of the present system was therefore 
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conducted; this analysis is reported in Section II. Sources of error which 
limit the accuracy of the readout system are identified; three sources 
wherein improvements are needed are sources which produce uncertainties 
in pattern line edge definition. 

With pattern line edge uncertainty identified as the primary limitation of 
ESVG readout accuracy, consideration is given to steps which can reduce 
the observed uncertainty. These steps fall into two categories: steps 
which reduce the effect or amplification of the error through the train of 
readout logic and the necessary calculations; and steps which reduce the 
uncertainty of the pattern line edge itself. 

Sections III and IV consider means of reducing the amplication of error. 

With a modest increase in the region of the rotor which is patterned and a 
change of readout strategy to one which utilizes three pickoffs when they 
view the pattern, a significant reduction of overall error amplification 
can be expected. Section V deals with prospective methods of defining 
the pattern edge more precisely. Some methods appear to be encouraging 
and should be pursued with laboratory experimentation. 

To accurately define the colatitude angle with the presently used approach, 
the longitude angle separating the points of the pattern lines viewed by the 
pickoff must be precisely defined. Section VI compares the advantages of 
two methods of measuring that angle -- an all-digital clock system which 
is the method presently employed and a hybrid phaselock system. If the 
phase discrimination can be made sufficiently sensitive, the hybrid system 
can be compatible with projected system requirements and offers an attrac- 
tive bonus -- the prospect of obtaining rate information from the ESVG. 

Section VII describes the effects of dynamic environments on the ESVG 
readout. With consideration for their effects in readout system logic design 
and computer compensation, the effects of static acceleration, vibration, 
and angular rotation can be reduced to tolerable levels. 
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SECTION II 

READOUT ERRORS OF EXISTING ESVG* S 


To accomplish an improvement in the readout accuracy of an ESVG, one 
must ascertain which of the individual error sources will significantly im- 
prove overall accuracy when the magnitude of these errors is reduced. Sub- 
sequent effort can then be concentrated on the sources so defined. An error 
analysis of existing gyro readout systems was therefore conducted to delin- 
eate the potentially fertile areas of improvement. 

The analysis begins by defining the physical sources at the fabrication, piece 
part, or subsystem level and describing its effect on readout accuracy. 
Methods of controlling and/or compensating these errors are discussed. 
Foreseeable improvements are indicated, and a breakdown of what the errors 
must be to achieve a one- second readout accuracy is provided. The areas 
where improvements are required or are of greatest benefit are then defined. 


ERRORS IN THE ROTOR PATTERN 


Pattern Application Error 

This is the low-frequency deviation of the application tool from the desired 
position of the pattern caused by an error in the master pattern (if a master 
is used) and by play in the tool driving linkage. Generally, a master is 
required for nonplanar patterns such as the cosine pattern. The deviation of 
the solid line from the null ordinate in Figure 1 illustrates the effect of this 
error when the nominal mechanization is used for calculating the direction 
cosine. 

This error is compensated in large part by altering the mechanization equa- 
tion. In the case of the cosine pattern, the nominal mechanization equation 
is 


cos 6 = A + B 0 
o o 

To compensate for application error, the altered equation could take the 
polynomial form 

n 

k=0 

with P q and P^ deviating slightly from A q and B q , respectively. All coeffi- 
cients are determined by a least- squares curve fit to calibration data. 


cos Q 


■1 
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Figure 1. Typical ESVG Rotor Optical Calibration 






Spin Axis - Pattern Axis Misalignment 

Associated with ESVG rotor patterns is an axis which should be aligned 
parallel to the spin axis. For several reasons, these axes may not be par- 
allel; consequently, a readout error can arise if coincidence is assumed. 

The error is a low-frequency error and can be compensated in large part by 

altering the nominal mechanization coefficients A and B . 

o o 

The misalignment of the pattern axis can be precisely determined in the pre- 
sent system by ascertaining that meridional angle which is simultaneously 
read on all three pickoff outputs (denoting that the spin axis angles with all 
three pickoff axes are equal). Geometrical considerations and knowledge of 
the location of the pickoff axes enable a calculation of the direction cosine 
associated with that angle. Two such angles can be determined, which are 
located symmetrically about the equator. Correlation of these results with 
the calibration data from which Figure 1 is plotted yields a calculation of the 
pattern axis misalignment. 


Line Edge Uncertainty 

A severe problem in any pattern application process is line edge control. In 
a masking technique, for instance, the mask which defines the line edge can 
be torn when the pattern line is defined on the surface or when the mask is 
stripped from the rotor. A grit-blast technique causes uncertainties because 
of the dispersion of the blast. Since the grit-blast method offers a distinct 
advantage of fewer steps, this method has been employed in patterning most 
rotors. Rotor noise, which is discussed later, produces an effect identical 
to edge dispersion; these contributions cannot be separated. The resulting 
error is a random function from rotor to rotor with characteristic frequencies 
ranging up to the 100th harmonic of the colatitude angle. 

Higher-frequency uncertainties can occur in the case of a nonplanar pattern 
cosine in which the application tool drive is servoed to a quantized master 
pattern line. This contributor produces smaller error than the lower-fre- 
quency contributors but can remain a significant factor by itself. The high- 
frequency characteristic of this error renders it impractical for compensa- 
tion; but, if readout occurs while the direction cosines are changing slightly, 
this error can be reduced to a certain extent by averaging. 

The scatter of the data points about the solid line in Figure 1 shows the effect 
of line edge uncertainties. Let xp denote the scatter observed in the meridi- 
onal angle measurement between pattern lines. From a statistical standpoint, 

an accumulated scatter of xp would be produced by i/z/Vi"" of scatter in each 
line. Since it is a measurement projected along a constant latitude line, it 
is an oblique projection of the line uncertainty (see Figure 2). The normal 
projection, e, can be calculated by 

e = -=£ sin 

Vi 


tan 


•1 d0 
d0 
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Figure 2. Amplification of Line Edge Uncertainty by 
Pattern Slope 


For the present gyro. 


de = 0. 522 
d0 sin 0 

Figure 2 shows how small slopes can amplify line edge errors in the longitude 
measurement. At a sacrifice in resolution afforded by small slopes, this 
amplification can be reduced by increasing the slope. 


ERRORS CONTRIBUTED BY THE ROTOR 


Reflectivity Variations 

The signal level received by the optical sensor is a function of the rotor re- 
flectivity. If triggering occurs at a fixed signal level, the fraction of the 
pickoff field of view that must be exposed to the pattern line would depend 
on the reflectivity of the unpatterned surface. The quality of the surface 
finish is a determining factor of reflectivity. Low-frequency variations can 
occur as a result of pqjishing variations on the sintered surface. 

Variations in reflectivity have not been measured on recently built gyros. It 
has been estimated that such variations are less than 10 percent of the nomi- 
nal signal. With the triggering voltage set at 50 percent of the nominal sig- 
nal, a 10-percent variation would cause up to 15 seconds of colatitude angle 
error. Compensation is provided by triggering at a point in proportion to 
the peak-to-peak signal, which reduces the error to that caused by nonlineari- 
ties in the sensor curve and the triggering electronics. 
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Rotor Noise 


Local variations in reflectivity occur as a result of surface perturbations and 
arcing between the rotor and the electrodes, giving rise to a noise component 
in the pickoff signal. Most of the noise pulses are electronically filtered by 
establishing a triggering threshold and pulse width discrimination. Very large 
pulses cannot be distinguished from pattern pulses and cause false triggering. 
Equally serious from an accuracy standpoint are marks and perturbations >on 
or adjacent to the pattern line. Their contribution of noise cannot be sepa- 
rated from the pattern signal and is not distinguished from pattern line edge 
uncertainties by the optics. 

An estimate of the error caused by rotor noise can be obtained by comparing 
the scatter in calibration data taken under gyro operating conditions with the 
scatter in the microscope calibration data. Such a comparison indicates that 
the rotor noise in one gyro contributed 22 seconds of error in the latitude 
angle. Like line edge uncertainties, this error source causes deviations 
about the solid line of Figure 1. 

Centrifugal Distortion 

Stress analyses of thin spherical shells'*' have predicted that the beryllium 
rotor radial contraction at the pole is approximately the same as the radial 
expansion at the equator as the rotor speed increases. One can look upon the 
strain as a migration of the rotor surface toward the equator; such a migra- 
tion causes an error in the readout. 

If the migration that occurs in the rotor can be reduced from a three-dimen- 
sion problem to a two- dimension problem, a quantitative calculation of the 
angular distortion can be calculated quite simply by utilizing the results of 
the stress analysis and tables of the elliptic functions and their integrals. 

For a two-inch rotor at a speed of 700 rps, the maximum migration predicted 
by this method is three seconds at colatitudes of 45° and 135°. 

Compensation can be implemented for centrifugal distortion by altering the 
mechanization equation appropriately. However, if the original calibration 
data is collected at operating speed and regressed for coefficients, the dis- 
tortion will automatically be taken into account. 

Nutation 

Readout errors can occur if the rotor nutation is not completely damped out 
during the starting process. In Honeywell ESVG history, nutation has not 
been a limiting factor in readout accuracy, either in the gimbaled or strap- 
down mode. If nutation does become a problem, the error can be reduced 
by averaging successive revolutions over the period of the nutation. 


"'■Refer to the Appendix of Volume II of this report. 


9 



OPTICAL SOURCES OF ERROR 


Illumination Variations 

If the field of view is not uniformly illuminated, the portion of the field occu- 
pied by the pattern line to cause triggering is a function of the pattern angle 
in the field. To reduce this effect, great care has been taken in designing 
and manufacturing pickoff lamps having uniform luminosity. Despite this 
effort, small variations have been known to remain. 

Field- of- View Irregularities 

If the field of view is not circular, the triggering point is dependent on field - 
of -view geometry and the angle of the pattern line in the field. At the assem- 
bly level, this source of error cannot be distinguished from illumination 
variations. Their combined effect can be observed by rotating the case of an 
operating gyro about the pickoff axis and noting changes in the phase angle as 
it rotates. 

Theoretically, a center of illumination can be defined, the location of which 
depends on the combined effect of illumination variations and field- of- view 
irregularities. If this axis is coincident with the desired location of the pick- 
off axis, readout error can be reduced to zero by triggering at 50 percent of 
the peak-to-peak signal. Other factors may render one or more of these con- 
ditions unattainable. If this is the case, compensation can be implemented 
by introducing a correction to the calculated direction cosine that is a func- 
tion of the pattern line angle of approach. 

The pickoff alignment procedure used on recently built gyros is one in which 
angle -of- approach sensitive errors can be a significant source of uncertainty 
in the alignment data. When effects which can be attributed to pickoff mis- 
alignment and spin axis drift are taken out of the data, a peak-to-peak varia- 
tion of 10 seconds in the colatitude angle has been noted. Since other sources 
of error contribute to the noise, it is difficult to discern how much of the noise 
is caused by illumination variations and field- of- view uncertainties. 

Lamp Luminosity Variations 

The luminosity of the lamp changes if the power supply fluctuates. Lamp 
efficiency can be expected to change as the lamp ages. Such changes produce 
changes in the signal level. If triggering occurs at a fixed level, an error in 
the readout angle from lamp power variations can occur. Power supply tran- 
sients can be controlled to some extent with regulation equipment. Lamp 
life is prolonged by extending the turn-on time and operating them at pru- 
dently selected color temperature levels. Extensive life tests have been con- 
ducted on pickoff lamps to determine how life changes with color temperature, 
to select the optimum operating temperature. 
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Electronic compensation is provided for any changes in power that do occur 
by the height-insensitive trigger in recent gyros. The triggering point is 
regulated to occur at a specified proportion of the peak-to-peak signal level. 


ELECTRONIC SOURCES OF ERROR 


Rise Time Variations 

Random variations in rise times of the signal-processing electronic compo- 
nents is a source of noise in the readout data. The vendor literature claims 
a rise time of 0. 2 microsecond for the photosensor used in one of the present 
gyros, and the rise time of the preamp is 0. 4 microsecond. The rise time 
of the sensor-preamp combination is taken to be the rss of the component 
rise times, or 0. 45 microsecond. The input pulse at a rotor speed of 700 
rps is 0. 9 microsecond, with a peak-to-peak variation of 10 percent occurring 
due to the change of pattern slope with latitude. An electronic variation of 
five percent would produce a 0. 01 microsecond error in the time interval 
measurement between pulses, which produces an error in the colatitude angle 
measurement of approximately three seconds. 

This error can be reduced by cutting down the rotor speed, averaging over 
several revolutions, and using higher-speed electronics components. Higher- 
speed components are available at a sacrifice in gain; however, one may take 
advantage of recent advancements to increase electronic speed without sacri- 
ficing gain and noise. It is believed that this error can be reduced to less 
than one second with present-day technology. In addition, the expected error 
from these variations can be reduced by averaging. 

Amplifier Noise 

Amplifier noise, like rotor noise, produces uncertainties in the readout mea- 
surement. A dominant source of noise is electromagnetic pickup of the sus- 
pension bridge frequencies by the electronics and the transmission lines. 

Triggering Level Error 

If triggering does not occur at 50 percent of the peak-to-peak signal, the 
pulse defining the crossing of a line edge occurs at some point other than the 
center of illumination. In addition to failing to null errors caused by 
illumination variations and field-of-view irregularities, a predictable pulse- 
timing bias that is a function of pattern slope occurs. It will be shown in 
Section III that the pattern can be designed to minimize the effects of this bias. 

Triggering Uncertainty 

In addition to the uncertainties already present in the input signal, intrinsic 
uncertainties can be developed by the pulse-shaping electronics of the height- 
insensitive trigger. For example, the output rise time of the present cir- 
cuit is of the order of 0. 1 microsecond. At the rotor speed for which it is 
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designed (200 rps), a gating error of 0. 1 microsecond would produce an 
error of about six seconds. At an operating speed of 700 rps, the same 
gating error causes an error of 21 seconds. 

Like the rise time uncertainties, triggering uncertainty can be reduced by 
the use of faster electronics, by averaging, and by operating at a lower rotor 
speed. Present electronic technologies would permit a higher- speed trigger 
circuit so that errors can be reduced to less than one second at a speed of 
200 rps. 


Phase Angle Measurement Error 

A digital method of phase angle measurement is employed in testing the pres- 
ent system. Clock pulses are counted, with the counter gated by the pulses 
from the readout system. In this system, three potential sources of error 
are considered. 

• Clock instability - No data is available concerning the 
clock stability. Test experience supports the belief 
that clock instabilities produce insignificant error com- 
pared to other sources. 

• Counter uncertainty - The quantized nature of a digital 
measurement is no more accurate than the least count 
of the counter. In the 100 me counters now used, the 
least count represents an error of 0. 01 microsecond 
per revolution, or two seconds of latitude angle error at 
700 rps. 

• Gating circuits - It is important that the gating circuits 
be sufficiently fast and uniform in their characteristics 
to be commensurate with the accuracy requirements of 
the gyro. The circuits presently used have been shown 
to introduce no more error than that introduced by the 
quantization error of the counter. 

All of the electronic sources of error are random and have high character- 
istic frequencies. The expected errors of all can be reduced by averaging 
over several revolutions. If the same approach is used in measuring phase 
angles in future systems, the use of the present electronic techniques, but 
with faster components and circuits, appears to be adequate for future ac- 
curacy requirements. 


OTHER SOURCES OF ERROR 


Pickoff Alignment 

To utilize the readout information, the location of the pickoff axis must be 
known to an accuracy commensurate with the system requirements. Pickoff 
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alignment errors are presently determined with pick off output data as the 
case of an operating gyro is rotated about that axis along which each pickoff 
axis is to coincide. The accuracy of the determination, therefore, depends 
on the accuracy of the readout system. 

In recently built gyros, the alignment of pickoffs can be adjusted so that the 
pickoff axis can nominally be aligned coincident with the desired axis. The 
accuracy of the determination depends on the range of the pattern used. The 
range decreases in proportion to the pickoff axis- desired axis angle; as the 
range becomes more restricted, the accuracy of the determination is improv- 
ed. The present accuracy limit is five seconds, the noise attributed to pick- 
off electronics and to illumination variations and field-of-view irregularities. 

Pickoff Alignment Stability 

The stability of the alignment generally is a function of temperature, accel- 
eration environment, and environment history. One can introduce compen- 
sation if a repeatable sensitivity to temperature is discovered. In other re- 
spects, the alignment stability depends on construction materials and mech- 
anical design of the pickoffs and their mounting. No data have been collected 
to evaluate the alignment stability or its sensitivity to temperature. 


Computational 


System computational accuracy depends on the accuracy of the direction 
cosines. In the present system, two direction cosines are determined dir- 
ectly from the pickoff signals and the pattern mechanization equation. A third 
cosine is calculated from pickoff alignment information and the identity which 
direction cosines of a unit vector satisfy in an orthogonal coordinate system: 

2 2 2 
cos 0^ + cos $2 + cos = 1 

The error generated in the third direction cosine from readout errors in the 
first two direction cosines can be estimated from the implicit differential 
of the identity. 


cos 0 1 + cos ©2 e 2 

cos 0o ~ 


where 


e i = A(cos 0^ = - sin 9 ^ A0.^ 

Not only is the error in the third direction cosine a function of the error in 
the other two cosines, but also of the cosines themselves. The range of the 
cosines from pickoff information is governed by the patterned region of the 
rotor. In the present system, 

I co s 0 I 
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From the orthogonality identity, it follows that the the magnitude of the third 
direction cosine can range from 0 to 1. Although it is ordinarily avoided in 
gyro operations, the possibility of a very small value for a third direction 
cosine and its effect on computational errors cannot be overlooked. 

The computational error here is a compounding of other readout errors that 
is reduced by reducing the other errors. The compounding effect can be re- 
duced by taking it into account in planning the operation of the readout system. 
This is considered in greater detail in Section IV. 


OBSERVATIONAL EFFECTS OF READOUT ERRORS 


The sources of error produce effects in the readout which follow distinctive 
patterns. The error sources may be classified according to the patterns 
observed in their effects. Often these patterns indicate how the effects of 
these errors may be compensated or reduced. 


Nonrandom Changes in Pickoff Signal and Signal-to-Noise Ratio 
Error sources which produce these effects are 

• Rotor reflectivity variations 

• Lamp lumininosity variations 

• Rotor noise isolated from the pattern 

• Illumination variations 

• Field- of- view irregularities 

The effects of the first two sources are nulled by the height-insensitive 
trigger (HIT). The operation of the HIT is discussed later. Most of the 
rotor noise is filtered out by establishing a triggering threshold and using 
a pulse width discriminator. The last two sources produce an irregular 
pulse shape that can be nulled in principle by aligning the center of illumi- 
nation to the desired readout axis and triggering at the 50-percent point of 
the pulse. 

Errors Which Correlate With Spin Axis-Pickoff Axis Angle 

Error sources whose effects can be correlated with pickoff axis colatitude 
are 

• Pattern application error 

• Pattern axis- spin axis misalignment 
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• Centrifugal distortion 

• Line edge uncertainties 

• Rotor noise at or adjacent to the pattern edge 

Compensation for these errors can be implemented by modifying the mech- 
anization equation for the direction cosine. Most of the error produced by 
the first three sources can be compensated by the addition of a few terms 
to the nominal mechanization. Many more terms must be added to handle 
line edge uncertainties and rotor noise effects. While these effects are ran- 
dom in the spatial domain, they are bounded; a potentially feasible method of 
compensation would be a Fourier expansion with the measured phase angle 
employed as the independent variable. The limit to which such an expansion 
can be carried would be defined by the resolution of the calibration data, the 
allotted computer capacity for the mechanization, or the point at which the 
benefits gained with the step wise addition of terms that correlate with the 
data no longer outweigh the inconvenience of carrying that additional term. 

Errors Sensitive to Pattern Line Angle of Approach 
Error sources whose effects are sensitive to angle of approach are 

• Illumination variations 

• Field-of-view irregularities 

• Pickoff misalignment 

Compensation for these effects can be introduced by subsequent mathematical 
operations on the direction cosines obtained from the mechanization equation. 
Pickoff misalignment is included here because the technique of compensation 
is similar. Moreover, the most significant effects of illumination variations 
and field-of-view irregularities cannot be separated from misalignment 
effects in an assembled gyro. 

Cyclic and Random Effects 
Error sources producing these effects are 

• Rotor nutation 

• Rise time variations 

• Amplifier noise 

• Triggering uncertainty 

• Phase -angle measurement error 
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The effect of rotor nutation can be nulled by averaging over the number of 
revolutions required for the spin vector to precess one revolution about the 
maximum moment axis. The expected error of the other effects is reduced 
by averaging (summing) over several revolutions of the rotor. 

Except for nutation, these error sources are electronic, the effect of which 
is a random error in the time axis. The effects of these errors can be 
reduced by other means as well as by averaging. The error itself can be 
reduced with the use of faster electronics. Such a reduction is within the 
scope of present technology. Another way to reduce it is to increase the mea- 
sured time intervals which can be accomplished by running the gyro at a 
lower speed. 


APPROACH TO ERROR COMPENSATION 

Four approaches to compensation are thus suggested as means to reduce the 
effects of errors without reducing the errors themselves. The historical 
development of these approaches is discussed in the order by which they 
would be employed. 


Height-Insensitive Trigger (HIT) 

Reflectivity variations and changes in lamp output appear as short-term and 
long-term changes in pulse height. If the trigger occurs at a fixed voltage, 
a time interval error (AT) results, as shown in Figure 3. 


SIGNAL 

VOLTAGE 


TRIGGER 

LEVEL 



Figure 3. Error From Fixed Trigger Level 
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Since the basic electronic trigger circuits are voltage- sensitive, some method 
of compensation is desired to remove this error source. One method is by 
compensating the trigger device to render it sensitive to a fixed proportion of 
height rather than a fixed voltage. This is the purpose of the HIT. The block 
diagram of the HIT is shown in Figure 4, and the waveforms are shown in 
Figure 5. The basic operation is as follows. 

The input pulse is amplified; the amplified pulse is sent to a peak level sen- 
sor that samples the peak voltage (flat top) of the input pulse and holds it. 

The amplified input pulse is then subtracted from one-half of the output of the 
peak level sensor. The zero-volt crossing of the resultant signal is a fixed 
percentage (nominally 50 percent) of the incoming pulse height. The zero- 
crossing signal is amplified to increase the slope of the zero crossing, and 
the zero crossing is then sensed by a tunnel diode threshold detector circuit. 
The trailing edge of the output pulse is the time mark used to control the 
time interval counters. The input pulse is differentiated and clipped, ampli- 
fied, and used to discharge the voltage level sensor so it can accurately 
sample the peak voltage of the next input pulse. 

Time Interval Measurement Averaging 

The meridional angle of the pattern is determined by calculating the ratio of 
the time interval measurement between the pattern lines and the time inter- 
val measurement of the rotor period. To obtain an average value of the 
meridional angle, both time interval measurements are summed over the de- 
sired number of revolutions. The logic necessary to perform the summing 
is built into digital data processing circuits that gate the high-frequency 
pulse counters used in measuring the time intervals and record the counter 
data on punched paper tape. A selector control enables an operator to aver- 
age over 1, 2, 4, 8, or 10 revolutions, the number chosen depending gen- 
erally on the nature of the test to be conducted and the rotor speed for the 
test. 


Compensation in the Direction Cosine Equation 

When a great circle pattern is used, the nominal mechanization equation 
takes the form: 

cos 6 = cos (cot" 1 (K cos 0) 


or 

COB e-- ssa - - 

Vl + cos (j) 

where 0 is the measured meridional angle, and K is a constant determined 
by the slope of the pattern line with respect to the spin axis. It is apparent 
that to mechanize this equation to the required accuracy in a vehicle-mounted 
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NOTE: CIRCLEO LETTERS REFER TO WAVEFORMS 
SHOWN IN FIGURE 5 

Figure 4. Height - Insensitive Trigger Block Diagram 
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Figure 5. Height - Insensitive Trigger Wave Forms 
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computer would require more arithmetic steps and memory bits than that 
required for the cosine pattern. The memory requirements of these patterns 
are compared in Section III. 

Computer compensation for cosine pattern application error and pattern axis- 
spin axis misalignment has been implemented by the addition of a few terms 
to the mechanization equation and regressing for all coefficients. Compen- 
sation for the higher-frequency errors (line edge uncertainties and rotor 
noise) has not yet been attempted. Experience in gathering calibration data 
from an operating gyro indicates that great care must be paid to procedural 
details and data gathering techniques in order to obtain calibration data in 
which uncertainties produced by spin vector drift does not exceed the accuracy 
requirements of the calibration. 

Once the data collection procedure is perfected, the problem of characterizing 
the higher-frequency errors remains. These effects are random (in the sense 
that they cannot be predicted a priori) but repeatable as a function of spin 
axis-pickoff axis angle. In addition, these errors are expected to be bounded; 
a potentially feasible method of compensation is a Fourier expansion which 
employs the meridional angle as the independent variable. While a Fourier 
form of this characterization is in principle not difficult to find, the imple- 
mentation of such a mechanization in a vehicle computer would require a con- 
siderable amount of memory. 

Compensation for Angle of Approach 

Compensation for pickoff misalignment has been implemented in a system of 
computer programs which reduces and analyzes drift data from the present 
strapdown ESVG. The location of the pickoffs in case axis coordinates is 
determined from the alignment calibration data. For each readout mode 
(combination of pickoffs read), an intermediate system is defined. One axis 

A 

of the coordinate system, X 1 , is coincident with one pickoff axis; a second 

A ■ 

axis, X 2 , is normal to the first but in the plane defined by the two pickoff 
axes. The angle between this axis and the second pickoff is equal to the non- 
orthogonality between the two pickoff axes. The third coordinate axis, Xg, 

is normal to the other two. Spin vector components in these coordinates are 
calculated as follows: 

S • = cos 0 1 


s- X 2 


cos eg ~ cos cos T 
sin t 


S 


X. 


= ±4 - (S 


V 


(S 


V 
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where 8. is the spin vector colatitude angle with the ith pickoff axis and t 

^ A A 

is the angle between the pickoff axes. The sign of S • is determined 

from continuity criteria. Spin vector components in case axis coordinates 
are then calculated from the transformation matrix of the intermediate set. 
Second and higher harmonics of the angle of approach have been noted in 
alignment calibration data collected from present gyros. The causes of 
these harmonics have not been fully investigated; in principle, there are three 
possibilities: field- of- view irregularities, illumination variations, and spin 
vector drift. Drift can be separated from the other causes by collecting 
calibration data with the gyro case rotating in each direction about the refer- 
ence axis and reducing the data from each simultaneously. 

If field-of-view irregularities and illumination variations are shown to pro- 
duce errors which correlate with higher harmonics of the angle of approach 
and which must be compensated, the method of compensation would likely 
take a different form from that described above for pickoff alignment. The 
direction of the spin vector projection on the plane normal to the pickoff axis 
would be determined from information provided by the other pickoffs. The 
angle of approach would be determined from this result and from the pattern 
slope at the observed colatitude. The correction, a function of the angle of 
approach, can then be calculated. 


CRITICAL AREAS OF REQUIRED IMPROVEMENT 


With the aid of Table I, the critical areas where improvements must be 
made can now be defined. Estimates of the error in present systems without 
compensation are given where data is not available. The data which is on 
hand is, however, classified; this data is furnished in Volume V of this 
report. Characteristic dependencies are listed, so that the error frequen- 
cies for a given system application may be calculated. Where compensation 
has been provided, the estimated error after compensation in the present 
system is indicated. 

Error estimates of an improved readout system are also listed. In all 
cases, the applicable form of compensation is assumed. The only source not 
considered capable of compensation is the master quantization error of a 
cosine pattern. Improvements in pattern, rotor, and optic sources of error 
are credited to the use of faster electronics which yield greater resolution 
in the calibration data on which compensation is based. No limits were 
assumed on the number of terms in the mechanization equation. System 
error estimates are 7. 5 seconds for a system employing a cosine pattern 
and 4. 5 seconds for on© using a great circle pattern. The difference between 
them is the absence of the master quantization error when using a great cir- 
cle pattern. 

An error budget which provides a system error of one second is presented in 
the last column of the table. The comparison of this error budget with the 
estimates of the improved system is the basis for defin ing those areas where 
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improvement is needed. The greatest improvement is needed in the reduc- 
tion of error from pattern edge uncertainties. Methods of reducing its effect 
are considered in Section IV, while methods of reducing the error itself are 
discussed in Section V. All contributing sources of edge undertainties -- 
which includes dispersion encountered in applying the pattern, arc marks, 
and surface finish imperfections -- are considered. 

Optical sources of error are presently not serious contributors. Electronic 
sources and centrifugal distortion are serious contributors at present; how- 
ever, they can be reduced to tolerable levels by using faster circuits and 
operating the gyro at a lower speed. Recent drift performance data from a 
700 rps unit and a 200 rps unit indicate that both have about the same perfor- 
mance potential. From a readout accuracy standpoint, the lower speed is 
desirable; from the performance data one may conclude that an intermediate 
operating speed would provide better drift performance. In either case, 
faster electronics and centrifugal distortion are not expected to produce 
errors which exceed the accuracy requirements of future ESVG units. 

The electronic equipment used to measure phase angles of present units is 
too bulky to accompany vehicle-mounted ESVG' s. Section VI considers a 
system which may prove to be feasible for use with such units. 

Improvement is needed in the control and compensation of pickoff alignments. 
Also, some assessment is needed of the alignment stability and how it is 
affected by environments. Laboratory effort is needed to provide improve- 
ments and knowledge in this area. Such an effort was not pursued in this 
study because it was beyond the scope of the study. 

Computational compounding of errors has been shown to be a problem worthy 
of further attention. This compounding is a major point in the study of read- 
out modes in Section IV. 

In this error analysis, no assessment was made of the environmental effects 
on readout. The effects of three environmental factors - steady acceleration, 
vibration, and vehicle angular motion - are considered in Section VII. 
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SECTION III 

PATTERN DESIGN STUDY 


An example was given in Section II which illustrates how the magnitude of 
one source of readout error is affected by the slope of the pattern line on the 
rotor. In this section we shall examine whereby several sources of error 
produce effects sensitive to pattern slope and other parameters such as rotor 
speed and the rotor colatitude observed by the pickoff axis; it will be seen to 
what extent error considerations have influenced the present pattern design 
concept. 

Error sensitivity equations in terms of selected pattern parameters are 
derived and compared for three pattern types --great circle, cosine, and 
colatitude. Optimum parameters are then selected which minimize the 
effects of the error sources. From the error amplification standpoint, the 
colatitude pattern is favored because its direction cosine error sensitivity is 
uniform and it has the smallest rms value. However, other factors can 
influence the choice of pattern type; two such factors are the ease of appli- 
cation and the mechanization equation computer memory size requirements. 
The great circle pattern is the easiest to apply to the rotor, while the cosine 
pattern requires the smallest memory size. 


PATTERN PARAMETERS 


The patterns considered here are defined by three parameters; pattern type, 
pattern range, and number of pattern lines. For any specific application, 
pattern range is likely to be determined by other considerations. However, a 
given pattern range may influence the choice of pattern type and the number 
of lines. 


Types of Pattern 

Patterns which have been applied to ESVG rotors to date are of two types: 
great circle and cosine. These two pattern types probably represent the 
extremes in difficulty of application to the rotor and in computer require- 
ments for mechanizing the direction cosine. The great circle pattern is an 
easily applied pattern because it is a planar pattern. The cosine pattern is 
desirable because it has a point-slope mechanization for the direction cosine. 

A third pattern type which may be desirable is one in which the colatitude 
angle is provided by a point-slope mechanization. This pattern shall be called 
the colatitude pattern. 

For a readout system in which a phase angle is measured and used to calculate 
the direction of the momentum vector, these three pattern types provide a 
reasonable diversification of features associated with many types of patterns 
compatible with this system. An error analysis of these features may then 
point out a desirable pattern or the desirable features of each pattern. 
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Pattern Range 


The pattern range is determined by whether or not all-attitude readout is 
required and the number and arrangement of the piekoffs. In most applica- 
tions, all-attitude readout capability is required; in certain applications, a 
limited-attitude capability may be sufficient, or all-attitude capability may be 
sacrificed for greater readout accuracy in the region of interest. For this 
study, it is worthwhile to consider the pattern range as a variable. 

Quantity of Pattern Lines 

In gyros which have been patterned to date, three lines have been employed — 
two pattern lines from which phase-angle information is derived and one logic 
line for separating pulse trains in the pickoff signal. The use of more than 
three lines has the advantage of reducing the effect of random errors by an 
averaging process. Other advantages may become apparent as the pattern 
study is pursued. 


NOMINAL PATTERN CHARACTERISTICS 


The direction cosine function has odd symmetry about the equator. Since good 
accuracy for all attitudes is of primary concern, there is no reason to con- 
sider error requirements which are not symmetric about the equator. The 
patterns to be considered in this study will, therefore, have the symmetry 
required to produce odd symmetry in the mechanization equation. The great 
circle and the cosine pattern which have been employed on gyros exhibit this 
symmetry, and the colatitude pattern to be considered in this study will also 
satisfy this symmetry criterion. 


The number of pattern lines is defined by the variable, n, which denotes the 
number of pairs of lines. At present, no consideration will be given to the 
logic line. A pair will consist of two adjacent lines, each of which covers the 
entire colatitude range of the patterned region of the rotor. Each line of a 
pair is a mirror image of the other with respect to a meridional line; i. e. , 
they have equal but opposite slopes at each colatitude point. It will be shown 
later that this condition minimizes the effect of a potentially significant source 
of error. 

Suppose that the lines of the first pair respectively satisfy the equations 

e = f u ( 0 ) 
e = f 12 ( 0 ) 
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The equational statement of this symmetry is 


d 0 d 0 

Each pair of pattern lines differs from another pair only in its meridional 
placement on the rotor. The pairs are equally distributed about the rotor in 
the longitudinal angle; no two lines overlap. Suppose that corresponding lines 
of the first two pairs, respectively, are segments of the following equations 

e = F n ( 0 ) 
e = f 21 ( 0 ) 

The criterion of equal distribution requires that 


F 21 <0> = F u (0 + -^) 

It is desirable to utilize as much meridional range for the pattern as possible. 
This desire, plus the criteria of slope symmetry and equal distribution of 
nonoverlapping lines, leads to the specification of two important mechaniza- 
tion parameters: 

• Each line is confined to a region having a meridional range 
of . The colatitude range of the line in the region will be 
defined by the lower bound and upper bound rr - 0 Q . 

• The equatorial intercepts of the pattern lines are equally 
distributed. The angle between intercepts of adjacent 

lines is, therefore, v = — . 

o n 


CALCULATION OF DIRECTION COSINE FROM PHASE-ANGLE DATA 


With the aid of the pair of pattern lines illustrated in Figure 6, the steps 
required to calculate the direction cosine may be outlined: 

• Determine the phase-angle v with the signal processing 
electronics 

• Calculate the meridional -angle p, from the pattern sym- 
metry characteristics 


Obtain the direction cosine from the equation of the 
pattern line 
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The details of the first step are treated in Section VI ; it is assumed here 
that v is known. The last step demands that the pattern equation be known 
explicitly; this will be derived later for each pattern type under considera- 
tion. The derivation of the angle n can be completed with the pattern charac- 
teristics that so far have been specified. 

From Figure 6, the angle (0 2 - 0^) can be expressed as a sum of three 

angles: 

i ji i 

v Q = 0 2 - 0i = (0 2 “ 02 ' + (02 " 0i ) + (0i " 0i) 



Figure 6. Pair of Pattern Lines on the Rotor 

Now, 

I 1 

V “ 0 2 — 0i 

which is known from phase-angle data. From the pattern line symmetry, 

^ 2 ” ^ 2 = 0.1 " 0.1 = ^ 
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Therefore 


If there are n pairs of pattern 

value of the set v, : 
k 


v 


Since 




This equation for p can then be substituted into the pattern equation to 
determine the colatitude or the direction cosine. 

The upper bound of p is useful for determining the parameters for the pattern 
slopes. Nominally, 

n 

0 ^ Z v, £ 2n 

k=l k 

The lower bound of this sum defines the upper bound of p. Therefore, 

p = JL 

max 2n 

The equational representations of the errors will define errors in terms of 
deviations, or expected deviations, of the meridional coordinate 0.' of each 

pattern line as sensed by the phase-angle measurement electronics. The 
sensitivity of the direction cosine to the errors is determined by tracing their 
effect through the chain of calculations which determine the direction cosine. 
Given that the pattern-mechanization equation can be written 

X = cos 0 = F (p) 

the sensitivity of X to an error in p is 

* x ■ § 


v Q = v + 2p 


ines, the phase angle provided is the average 


i n 

= l T 

n - i ~‘. 


v, 


n k“?i ' k 


V = 11 
o n 


= (tt “ X. v k> 


(i) 


t=i 
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The standard deviations of Ak and Ap are similarly related (with precautions 
taken to assure that standard deviations do not become negative): 


b (AM 


dF 

dp 


a (Ap) 


The sensitivity of p to an error in 0.’ can be similarly related. Let the 

standard deviation of the meridional-angle error be denoted by 0 (A 0.). Since 

all of the pattern lines are alike or symmetric, these standard deviations are 
the same. Since two pattern lines define each angle and since there are 

n pairs of lines. 


Thus, 


0 (A 


n 

£ 




= yin a(4 0 j) 


a (Ap) 


_1_ 

2n 


n 

a ( I] 
k=l 



a (A 0) 

y2n 


(2) 


MECHANIZATION EQUATIONS FOR THE PATTERN TYPES 


Great Circle Pattern 

The general equation for a great circle on a spherical surface is 

a sin 0 cos 0 + b sin 0 sin 0 + c cos 0 = 0 

where 0 is the colatitude angle; 0 is the longitudinal angle of the spherical 
polar coordinates; and a, b, c are the direction cosines of the axis normal 
to the plane of the great circle with respect to the X, Y, Z cartesian co- 
ordinates, respectively. 

The great circle pattern lines which would be applied to an ESVG rotor are 
segments of lines satisfying this equation. The parameter of these segments 
can be determined by combinations of a, b, c and the coordinates of the end 
points of the segments. 

For this design study, only one pattern line needs to be considered (which 
corresponds to Line 1 of Figure 6). The direction cosine solution from this 
great circle segment can be written 

a 1 cos p + b 1 sin p 

X. = cos 0 = - = 

V 2 o 

c.i + (a^ cos p + b^ sin p) 
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From the fact that p 


0 at 6 = , 


a l 


0 


Let be the slope of the pattern line at the equator: 


K 


1 




The direction cosine mechanization equation is, therefore. 


X 


Kj sin p 




2 . 2 
sin jx 


(3) 


The parameter can be more tangibly defined in terms of the pattern 
bounds with the following equivalent form: 

cot 6 = sinp (4) 


At the lower bound of the pattern region, 

TT 


^ = 2n 

e = e 

o 


Therefore, 


cot 0 


K 1 = 


sin 


TT 

2n 


(5) 


For the error analysis, the equation for the direction cosine derivative is 
needed: 


dX dX cos p 

d0 dp [1 + Kj^ sin 3 p] 3 ^ 3 

Equation (4) can be used for writing this derivative in terms of the colatitude 
angle. 
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= sin 3 0 Vk l 2 - cot 2 0 (6) 

From Equation (4), it is seen that K 1 and - K 1 define the upper and lower 

bounds of cot 0; thus, ^ is real for all values of 0 intercepted by the great 

circle. Combining Equations (2), (5), and (6) provides the equation relating 
the error sensitivity of the direction cosine to a (A 0): 

sin 3 0 / cot 2 0 P 

a(£\) = — - / 9 “ cot 9 a ( A 0> ( 7 > 

Vsin 2 ^ 

Cosine Pattern 

The general equation of the cosine pattern line is 

cos 0 = A 0 + B 

As it was developed for the great circle pattern, it is desirable to formulate 
parameters for the pattern characteristics by the equatorial intercept and its 
slope at the intercept. Again, only one line is considered: 

X = cos 0 = A^ p + 


for 


M- = 0, ©*«§■; 

therefore. 


The mechanization simplicity is demonstrated by the derivative of the cosine 


dX _ dX _ . 

d0 dp A 1 


( 8 ) 


As can be expected, this derivative is independent of the colatitude angle. 
The parameter A^ is determined by the coordinates of the pattern line at the 

colatitude lower bound of the pattern region: 

2n cos 0 

/V _ / r\\ 
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The pattern mechanization equation is, therefore. 


2n cos 0 

cos 0 = — — — 2 — fx (10) 

Combining (2), (8), and (9) yields the error sensitivity to o (A 0) for the 
cosine pattern: -, s 

~\f 2n cos 6 

o(A\) = — -2— o (A 0) (11) 

Colatitude Pattern 

The general equation of the colatitude pattern line is 

6 = M 0 + N (12) 

Again, following the same procedure, the equation for the line of interest is 

6 = p + 

For 

H = 0, 6 = ^ 


therefore. 


N = — 
iN l 2 


Therefore, 


X. = cos 6 = cos (M. p + - 7 j-) = - sin p 


(13) 


is the general pattern mechanization equation. Using the value of p at the 
pattern lower bound (p = ^ at 6 = 6 ) on the general form Equation (12) 
produces the result 


M x = 2n| 


6 


TT 


It should be noted that since 6 q < ~ is always negative 
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The derivative of the direction cosine function is 


■fr - - M i cos M i >* 

With Equation (13) and a series of trigonometric manipulations, one can show 
that 


^ = -M, sin 0 

dp 1 

Let 


The error sensitivity equation can now be written. 

a (AX.) = T sin 0 ~\/2n a(A 0) 


(14) 


(15) 


ERROR MINIMIZATION CRITERIA 


There are a number of error sources in which the system uncertainty from 
these errors is sensitive to certain rotor pattern characteristics. Four 
such characteristics are discussed and derived here. 


Sensitivity to a Phase-Angle Measurement Error 

Since the phase angle is the angle which is calculated directly from the pick- 
off signals, it is desirable to ascertain the sensitivity of the readout variables 
(direction cosines, or direction angles) to a phase-angle error. A generalized 
error equation can be derived in function notation. 

Let the direction cosine be the variable of interest for the moment. Given 
that 


e * Fj (0) 

the inverse function can be written 

0 = F 1 “ 1 ( 6 ) 


Now let 


X = cos 0 = cos [Fj (0)] = (0) 


34 



An error in the phase angle produces an error defined by 


d GA0) 

AX * — 1 A 0 = G / 1 (0) A 0 


It is desirable to state the error equation as a function of the colatitude angle 

0 . 


A X = G\ CF 1 -1 (0)] A 0 

An expected value for the magnitude of A 0 is the standard deviation of A 0, 
which shall be denoted by a(A 0). If A 0 is replaced by or (A0), an expectation 
of the magnitude of AX is the result. This entity shall be similarly expressed 
as 

a (AX) = \G' [F" 1 (0) ] | a (A0) (16) 

To assure a positive result for a (AX), the absolute value of the function is 
taken. In the case of a phase-angle measurement error, from electronic 
sources, a (A0) is independent of 0. 

Sensitivity to Line Edge Error 

Line edge error is a function of the patterning process; as such, the line edge 
deviation is measured normal to the line itself. Its projection in any other 
direction is amplified; therefore, the projection of line edge deviation along 
a parallel of latitude is a function of the slope of the line at that point 
(Figure 7). 



Figure 7. Geometry of Line Edge Error 
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The general equation for line edge error sensitivity can also be derived in 
terms of function notation. 

Referring to Figure 7, let e denote the expected value of the normal pro- 
jection of the line edge dispersion and d the resulting longitudinal dispersion, 
so that 


a (d) = e esc ^ 

where 

tan '^ = 7TJ- (17) 

Since 


esc 



1 

tan"* x]/ 


the longitudinal dispersion can be expressed in terms of the derivative 


0 < d > ■ e O 

d 0 



2 


The longitudinal error is related to a as follows: 


d = r sin 8 A 0 


so that 


cr (A0) = 


r sin 8 


r sin0 -T- r 
d 0 


1 1 + 


d 8 

id 0, 


(18) 


Later, it will be convenient to have a (A0) expressed in terms of the direc- 
tion cosine derivative: 


a (A 0) 



(19) 
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To determine the error in direction cosine, this equation is substituted into 
the equation for the phase angle error: 


A X 


r 1 


G, [F 


-1 


Jl + 


d e 
11 


sin 6 


d e 

d 0 


Now the derivative can be denoted by 


44 = F ' < 0 ) 

d 0 l 


so that 


sin 6 


d e 

d 0 


f* { - X) - - G i w > 


The final generalized equation for the direction cosine error is 


A x s f yi + CF; [ F^ 1 (0) 3 } 2 (20) 

(Since the absolute value of AX is of primary interest, signs have been dis- 
regarded,. ) 

Sensitivity to Triggering Level Error 

The pickoff signal is used to produce pulses for the phase-angle measurement 
electronics. This signal is changing as long as the pattern edge traverses 
the pickoff field of view, which represents an angle of the order of several 
minutes. Consequently, the triggering level of the signal must be pre- 
cisely defined. The ideal triggering point is the instant at which the line 
crosses the center of the field of view. 

It is worth-while to evaluate the dependence of errors caused by triggering 
level deviations on the pattern parameters. Let A denote the unpatterned 
surface area viewed by the pickoff (Figure 8). 
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a 2,2 2 . v 

A = tt p - (p ct - p sm a cos CT ) 

The signal output is proportional to the complement of this area. An expres- 
sion for the norm alized complement of this area is 




Figure 8. Pattern Edge in Field of View 

Let L be defined as the time at which the edge of the ith pattern line crosses 
the center of the field of view. Then 

h = ^ ,p - * - T i * G “ V * T i (22) 

i 

where T^ is the time required for the line edge to travel from the ege to the 
center of the field. If v is the linear speed of the line edge, 

X = - = - /oo\ 

i v cos j/a r sin 0 cos i/z. ' ' 
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Let the triggering point be defined as t. ' the instant when the normalized 
signal S is equal to ^ . Therefore, 



If the region of interest is restricted to 

0-4 ^ ri ^ 0.6 

a good approximation to Equation (24) is 


t.' - t. \ 

- V ~ - (25) 


The time interval between pulses generated by two pattern lines (call them 
lines 1 and 2) would be 

At (t|) = t 2 ' - t 2 ' =f (0. 5 - v) (T 2 - T x ) + t 2 - 


rr (0. 5 - x\) - 2 


If the trigger level is perfectly set at 

- 0, 5 

the desired time interval would be obtained; i. e. , 

At(0. 5) = t 2 - t 1 

The error in the time interval caused by an error in the trigger level setting 
is, therefore. 


e (At) = -y (0.5 - ^ (Tg - T^) 


in view of Equation (23). 


tt ( 0 . 5 - rpp / 1 

2 w r sin 6 cos \I/ 0 
s , & 


1 1 

cos 

J 


(26) 
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Now, if pattern lines with complementary slope functions [i. e, , i/^©) = 

-i//^(0)] are used, then cos i/Zg = cos and triggering level errors are 

eliminated. The correct time interval between pulses would be generated 
regardless of where the triggering level is set. Some distinct advantages 
develop as a result of eliminating this error. 

e Triggering levels may be adjusted to eliminate false 
triggering from arc marks, if necessary. 

• No compensation is necessary for the possibility that 
levels may vary among the pickoff s. Measures 
such as a separate mechanization equation for each 
pickoff are eliminated. 

For these reasons, pattern designs which are considered in this study will 
exhibit this slope symmetry between pattern line pairs. 


Sensitivity to Signal Rise Time 

The accuracy of the electronics which produces a pulse at the trigger level 
depends on the signal rate of rise at the triggering point. Differentiating the 
equation for the normalized pickoff signal Equation (21) with respect to time 
yields 


dS = dS dh = _ 2 / |_h_l 2 dh 
dt “ dh dt ~ ” np j p f dt 

From Equations (22) and (23) 


2 _ 

rtp 


sin a 


dh 

dt 


dh 

dt 



to r sin0 cos \p. 
s r i 


Therefore, 



2 . 

— sin q u r sm 0 cos \1/. 

TTp S r l 


(27) 


For the pattern analysis, only the pattern-sensitive parameters need to be 
carried explicitly - namely, 0 and \ p. Assuming that the readout error is in- 
versely proportional to the rate of rise, the longitudinal error generated by 
rate of rise attributed to pattern factors is proportional to e, where 

1 

' sin 0 cos 
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(It is no longer necessary to carry the subscript to for this analysis. ) 


By using Equation (17) and trigonometric relationships. 


Consequently, 


cos i// 


d 6 



e 


sin 6 


d 0 
d a 



A comparison of this Equation with Equation (18), in view of the fact that 
A 0 « e, shows that the cosine error function from rate of rise takes the 
same form as that from line edge dispersion. Since a pattern which reduces 
one of these errors reduces the other in the same way. Equation (19) will be 
considered to represent both factors in this study. 


Sensitivity of Readout Variables to Error Sources 

There are, therefore, two error-sensitivity equations from which pattern 
characteristics can be evaluated. The direction cosine error has already 
been derived. It is also desirable to know the colatitude angle error sensi- 
tivity to these effects as well as some idea as to how sensitive the third dir- 
ection cosine is to these errors. Both of these can be obtained from the 
basic direction cosine sensitivities. Since 

A X. = A (cos 6) = - sin 6 A 6 
the colatitude sensitivity (disregarding signs) is 

a (A 6) = — i — 5 - a (A X) = esc 6 a ( AX) (28) 

sin u 


From the orthogonality identity, 

X l 2 + S 2 + S 2 = 1 

therefore, 

X i A X i + X £ A X £ + X3 A X g = 0 


(29) 
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It is convenient in showing the third direction cosine sensitivity that only two 
direction cosines be considered. To accomplish this, it shall be assumed 
that ^2 zero > 

Then, 

AX 3 " XT AX i 


and 


therefore. 



AX 3 



AX 


1 


By putting = X and disregarding the signs, we find that 


a (AXg) = cot 9 a (AX) (30) 

Once the direction cosine sensitivity is derived, the colatitude sensitivity and 
third direction cosine sensitivities can be expressed directly through Equa- 
tions (28) and (30). 


CALCULATION OF ERROR SENSITIVITIES 


The sensitivity of the direction cosine equations for each pattern type can 
now be evaluated for two types of error. To evaluate the sensitivity to a 
phase-angle measurement error, let the standard deviation of the meridional 
angle error be given by 


cr(A0) = ut 

where t is the intrinsic timing error of the phase-angle electronics, and to 
is the rotor speed. To determine the sensitivity to line edge dispersion. 
Equation (19) and the direction cosine derivative for the pattern type under 
consideration are used. Each of these expressions for ct(A 0) is substituted 
into the error sensitivity equation derived for each pattern type. 

To determine optimum values of the number of pattern line pairs for given 
pattern ranges, an expression for cr(AX), which represents the error for 
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the entire pattern range, is desired, rather than the error at a given colati- 
tude. A convenient representation is the rms value of cr(AX) over the nominal 
range: 


CT = 


a 


rms 



1/2 


Great Circle Pattern 

For the great circle pattern, the sensitivity of the direction cosine to a tim- 
ing error is 


a (AX) = M- sin 3 9 


cot 3 9 


Viri 


sm 


_2 rr_ 
2 n 


;?■ - cot 3 6 


The rms value of the expected timing error is 


UT 


. 2 TT 

Sin 2H 


c o cot 9 cos 4 9 
iJL a J. O o 

*16 cot d o + TT - 26 


i 5 3 . 2 n 

|-g+ T sm 8 o 


2 2 
sin 0 O cos 9 C 


1 sin 9 cos 0 I 1 sin 3 0 sin 2 0 cos 2 0 ] 
• ° -2. LL o + 0 oJ 


IS rr-20 


8 


1/2 


The expected line edge dispersion equation is written in terms of a direction 
cosine derivative; the equation for the derivative is substituted before the 
final equation is obtained. Using Equations (5) and (6) in Equation (19), one 
obtains 


cr(A0) 


JL 

r 


— \/csc 9 + 


6 2 
sin 9 I cot 9 


sin 


2 rr 
2n 


cot 3 0 


(31) 
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The rms over the pattern range is 


a 


els 


b -qs { 


sin 


2 TT 


2n L 


2 

o 9 cos 0 

— cot 9 + 9fl — a 

4 o tt-20 


ircot0 + sin0 cos 0 
2 o o o 


(32) 


4 


7 - 


sin 2 0 cos 2 0 
o o_ 

TT - 20 


1/2 


A graph of each rms function is shown in Figures 9 and 10 for various values 
of 0 . The normalized timing error curves (Figure 9) show that the number 

of pattern line pairs should be kept to a minimum, while the normalized line 
edge error curves establish optimum values for n. In the ESVG, pattern 
line edge errors are expected to be greater than timing errors; the choice of 
n is, therefore, likely to be dictated by the optimum values for pattern line 
edge errors. The timing error contribution motivates a choice of one of the 
lower values of n which optimizes the pattern line edge errors. 


Cosine Pattern 


For the cosine pattern, the expected deviation from a phase-angle measure- 
ment is 


V 2n cos 0 

a (AX) = gj t - 


(33) 


Since Equation (33) is independent of the colatitude, it is equal to the rms 
value of that error. 


From Equations (8), (9), and (19), the expected error from the line edge 
dispersion is obtained. 


o (A0) 






V 


CSC 0 + 


TT 



4n cos 4 0 


Substitution of this expression into Equation (11) yields 


a(AX) 



2n cos 2 0 Q esc 2 0 


TT 



(34) 
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NUMBER OF PAIRS OF PATTERN LINES 


Various Pattern Ranges 



The normalized rms representation of this error is 


CT 

e 


_r 

e 


1 

2n 


+ 


“3 

4n eos 9 

o 

rr 2 sin 9 

o 


1/2 


( 35 ) 


Figures 11 and 12 are graphs of Equations (33) and (35), respectively, with 
the lower bound of the pattern range as a parameter. Again, the line edge 

deviation error shows that optimum choice of n can be made to minimize it, 
while the timing error motivates the use of a lower value of the optimums. 

Colatitude Pattern 


For the colatitude pattern, the expected phase-angle timing error is 


CT (AX) = w t T sin 9 2n 


(36) 


The rms value of this error over the pattern range is 

cr / ~ i sin 2 9 


CJT 


r-p 2 

T n 


1 + 


o 


TT - 2 9 


(37) 


The equation for the expected line edge dispersion is obtained from (14) and 
(19). 


ct( A0) = 


€ 

r sin 9 ~\J 


1 + 


1 


4 n 2 T 2 


Substitution into (15) yields 

o(AK) = 

Since this expression is independent of the colatitude, it also represents the 
rms of the error. Therefore, 


* e -f = V 2nT2 + ^ (38) 

Figures 13 and 14 are graphs of the error Equations (37) and (38), They show 
the same characteristics as the corresponding graphs for the great circle 
and cosine patterns. 
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COMPARISON OF PATTERN TYPES 


The three types of patterns which have been described can be compared to 
determine their relative utility in the ESVG. Of immediate interest are 
comparisons in three respects: amplification of timing and line edge disper- 
sion errors, ease of applying the pattern to the rotor, and ease of mechanizing 
each pattern equation in a computer. 


Comparison of Error Sensitivities 

A comparison of the graphs shows that the readout error introduced by the 
timing uncertainties is not a strong function of the type of pattern used. Ex- 
cept for the larger- ranged patterns and for small numbers of pattern lines, 
the curves of Figures 9, 11, and 13 coincide. The colatitude and cosine 
patterns tend to result in smaller errors from timing uncertainties than the 
great circle; however, the difference is not very great. 

The type of pattern has a greater influence on the expected error from line 
edge uncertainties, especially for larger-ranged patterns. The colatitude 
pattern produces the least expected ms error, while the cosine pattern pro- 
duces the greatest. A better comparison of these differences is obtained by 
selecting a common pattern range and a common optimum value for n and 
examining the expected error from pattern edge dispersion as a function of 
the colatitude angle. 

Since previous designs utilize a range of 90° (0 Q - 45°), it seems advisable 

to examine this range as a basis of comparison. The optimum value of n 
for all three pattern types is two. The normalized expected direction cosine 
error from line edge dispersion given by Equations (31), (34), and (38) is 
compared in Figure 15 for = 45° and n = 2. The normalized expected 

third direction cosine errors from these patterns, which are calculated using 
Equation (30), are compared in Figure 16. That the error vanishes at 
6 = 90° is not surprising when it is realized that under the assumptions 
leading to Equation (30) the magnitude of the third direction cosine is unity. 

From Figure 15, it is observed that the cosine pattern produces the least sen- 
sitivity to edge dispersion in the central pattern region, whereas the great 
circle produces the least at the pattern extremities. The colatitude pattern 
provides uniform sensitivity. If third direction cosine errors were not sen- 
sitive to the direction cosines themselves, the colatitude pattern would be 
the optimum pattern with respect to edge dispersion. 

With the presently used mode of obtaining three direction cosines from infor- 
mation given by two pickoffs, a pattern type that provides relatively small 
error amplification at the pattern extremities (the regions most likely read 
when the third direction cosine sensitivity is high) is the type which would 
likely provide the most uniform error amplification, if not the least amplifi- 
cation. However, other modes are considered in the next section; it will be 
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Figure 15. Comparison of Direction Cosine Errors from Pattern 
Edge Dispersion for Three Pattern Types 
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seen that the mode itself influences the amplification of edge errors. The 
optimum combination of pattern parameters and readout modes are better 
discussed there. 


Comparison of Ease of Application 

An important phase of the application process is that of defining the exact 
location of the line edge on the rotor. Since the great circle pattern is a 
planar pattern, its most distinct advantage is in the simplicity of defining • 
its line edge. Figure 17 illustrates how a great circle line segment can be 
defined. The rotor is rotated about the axis normal to the plane defining the 
great circle. A reference point, normal to this axis which intersects the 
center of the rotor, defines the great circle of which the pattern segment is 
a part. 

That this method of defining great circles is easily mechanized with hardware 
is the reason that great circle patterns were applied to the rotor when the 
concept of phase-angle measurement of pattern pulses was first utilized. The 
reference point is rigidly mounted on the apparatus; this rigid mounting eli- 
minates error sources that occur in defining the other pattern types. 

The definition of the other types of pattern lines is more difficult, since they 
are nonplanar. Figure 18 shows the apparatus presently used to apply a 
cosine pattern to the rotor. At the right is a double -size sphere on which a 
master pattern is scribed. The product rotor is mounted at the left. A pick- 
off views the master sphere; its mounting is servoed so that the pickoff is 
focused on the piaster pattern line. The angle of the mount is mechanically 
transferred to the grit nozzle mount at the left. The nozzle edge is the refer- 
ence point that defines the cosine line. This arrangement can be used for any 
pattern which is scribed on the master sphere. Errors in addition to those 
obtained by the great circle arrangement are errors existing in the master 
line and nonlinearities of the transfer linkage. 

A precise method of defining the pattern line on the master sphere is with a 
series of mathematically exact points. Errors in spindle scales can produce 
macroscopic deviations. Another problem is the quantization introduced with 
the use of points. The quantization is reduced by reducing the distance 
between points and averaging with a larger pickoff field of view. 

Other methods of defining the colatitude and cosine patterns can be devised by 
considering the mathematical characteristics of each. A method for defining 
a colatitude line is shown in Figure 19. The reference point is mounted on an 
arm which is rotated about an axis that lies in the rotor equatorial plane. The 
rotor is rotated simultaneously, such that the ratio of the two rotation angles 
is constant. Figure 20 is a diagram of a cosine line defined on the rotor. 

The reference point is translated in the direction parallel to rotor spin axis. 
The reference point displacement (z), in this direction from the equatorial 
plane, is proportional to the cosine of the colatitude coordinate of the plane 
parallel to the equatorial plane containing the reference point. If the ratio of 
displacement change to change in rotor rotation about the spin axis is kept 
constant, the locus of points at which a line from the reference point inter- 
secting the spin axis at right angles strikes the rotor surface defines a cosine 
line. 
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Figure 19. A Colatitude Pattern Applicator 




These methods of defining pattern lines to the desired accuracy with hardware 
is not necessarily within the current state of art. The arrangement indicated 
for defining the colatitude line does not appear to be too difficult. It is more 
practical than the arrangement for the cosine line definition. For this line, 
it is difficult to make the required transfer from the reference point-through 
a changing distance, t , to the surface whose normal to this line is also 
changing-with the required accuracy. This aspect renders this method very 
difficult to implement without first experimenting with new techniques. 

For ease of application, the great circle pattern is definitely the easiest to 
pattern to the required accuracy. The colatitude line is more difficult to 
apply to the required accuracy, but appears to be within reach. The cosine 
line is the most difficult of the three. Some laboratory research into the 
engineering feasibility of defining the cosine and colatitude lines with the 
methods discussed here or with similar methods is warranted. 


Comparison of Computer Mechanization Requirements 

To determine the computer mechanization requirements, each of the mech- 
anization equations must be broken down into the basic operations required 
for automatic computation. The basic mechanization equations are 

K. sin (jl 

k = 1 9 — (great circle) 

-yi + Kj sin 

k = A^p 
k =. -sin M^p 

where 


n 

' 2k <T1 ' Z 

k=l 

In addition to the equations, the accuracy requirements of the direction 
cosine and some estimate of the ranges of the variables and magnitudes of 
the constants are needed. 

To provide an angular accuracy of one second at the pattern range extremities, 
the direction cosine should be known accurately to the sixth decimal place. 

A binary representation consisting of 24 bits provides adequate accuracy and 
safety margin for computational error buildup. For this comparison, it is 
assumed that temporary storage locations and numerical constants contain 
24 binary bits. 


(cosine) 

(colatitude) 


v k) = 


TT 

2n 


v_ 
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The sequence of commands needed to complete the computation depends in 
part on how information is arranged on the accumulators and registers to 
carry out the arithmetic operations. It is convenient to choose a specific 
computer to provide these guidelines. The Packard-Bell 250 computer was 
chosen for this purpose. 

The computation of p_ is the first step and is independent of the pattern type. 

If one assumes that u is furnished the computer by the phase-angle measure- 
ment electronics, p is calculated by the steps shown in Table II. Although 
division by two is simpler by shifting the contents of the register one bit to 
the right, a division command was provided so that the input may be properly 
scaled. (This could be accomplished with a multiplication operation also. ) 


TABLE 1I.‘ COMMAND SEQUENCE FOR CALCULATING p 


OP Code 

— 

Location 

Nature of Operation 

CIB 


Clear Input Buffer 

TES 


Waits for signal indicating Buffer is 
loaded with v 

LAI 


Loads register A with V 

LDC 

C 1 

Loads register C with divisor 

DIV 


Quotient is formed in Register B 

STB 

X 1 

Stores B temporarily in 

LDA 

C 2 

Loads A with constant 

SUB 

X 1 

Subtracts contents of X^ from A 


>{c 

C. denotes a constant in memory 
Xj denotes a temporary storage location 

It is convenient to consider the pattern mechanization equations in increasing 
order of complexity. The simplest one is the cosine mechanization. To 
complete the calculation for X , the additional commands of Table III are 
needed. 

The computer may then proceed to other computations necessary for vehicle 
navigation. With the cosine pattern, a total of 12 commands, three constants, 
and one temporary storage location are needed in the computer memory to 
calculate the direction cosine. 
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TABLE HI. - COMMAND SEQUENCE FOR NOMINAL 
COSINE MECHANIZATION 


OP Code 

Location 

Nature of Operation 

IA.C 


Interchanges contents of A and C 
Registers 

LDB 

C 3 

Loads B with Multiplier A^ 

MUP 


Product is formed in A 

STA 

X 

Stores direction cosine 


To compensate for low-frequency deviations of the pattern line from the 
nominal location, the mechanization equation might take the form of a poly- 
nomial: 


X = P q + + P 2 |a 2 + + P m |a m 

It is of passing interest to ascertain the memory required for such a mechani- 
zation. A sequence of commands for this equation is given in Table IV. The 
bracketed steps would be repeated for each term in the polynomial beyond the 
second (linear) term. The total memory required, as a function of m , may 
be written as follows: 

Commands: 14 + 7 (m - 2) 

Temporary Storage: 2 

Constants: m + 2 

The colatitude mechanization is considered next as it is less complex than the 
great circle mechanization. The series expansion of the sine function is used 
to calculate X . To determine how many terms are needed to provide the 
desired accuracy, an estimate of the range of the product M^p is needed. 

Since the cosine of the colatitude is equal to the sine of the latitude, it is 
easily seen that the term (-M^) ig; equal to the latitude of the pickoff in 

radians. If 50° is a practical latitude limit on the pattern range, then 


|M 1 p I < 1 
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TABLE IV. - COMMAND SEQUENCE FOR POLYNOMIAL 
MECHANIZATION DIRECTION COSINE 


OP Code* 

Location 

Nature of Operation 

| fl 

IAC 


Places p in C 

STC 

X 1 

Stores p 

LDB 

C 4 

Loads B with P^ 

MUP 


Product P^ p formed in A 

ADD 

C 3 

Pq added to P^ p 


STA 

X 2 

Partial sum stored 


LDB 

X 1 

Loads B with p 


MUP. 


k “* 1 k 

p • P = P formed in A 


IAC 


k-1 k 

Replaces p with p in C 


LDB 

C 

k+3 

Loads B with P^ 


MUP 


Product P^ pk formed 


ADD 


Product added to partial sum 

STA 

X 

Stores direction cosine 


*The bracketed terms are repeated for each additional 
term in the polynomial. 
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The term containing gy- is the last term needed in the expansion to give the 

desired accuracy; therefore, five terms are needed. A sequence of com- 
mands providing this calculation is given in Table V. 

The bracketed terms are repeated four times to generate the sine function 
from the first five terms of the series. The memory required for the colati- 
tude pattern mechanization (for m = 5) is 

Commands: 16+7 (m-1) = 44 

Temporary Storage: 2 

Constants: m + 2 = 7 

The great circle utilizes a sine function also; however, the possible range of 
the argument is different: 


In!* 


TT 

2n 


To allow for the possibility that n = 1 , a possible value of for p is 

assumed. To obtain the desired accuracy, seven terms of the series for the 
sine function are required. The total number of commands required to form 
sin p can be determined by examining the command sequence for the colati- 
tude mechanization. Since p is not multiplied by a constant, the first three 
commands in Table V are not needed. The memory required to compute 
sin p (for m = 7) is: 

Commands: 13 + 7 (m-1) = 55 


Temporary Storage: 2 

Constants: m + 1 


8 


To complete the calculation, a square root must be formed. Generally, an 
iterative method is used. If the first estimate of the root is sufficiently pre- 
cise, Newton’ s process converges rapidly. To determine an equation for 
the first estimate and the number of iterations, an estimate of the magnitude 
of (K^ sin p) is needed. An examination of Equations (16) and (20) leads to 

the conclusion that 


I K. sin p 1 ^ cot 9 
l ' ' ” o 

If a practical range limit of 40° is assumed, the maximum value of this 
factor is 1.2. Thus, 

1 si + sin^p < 2. 5 
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TABLE V. - INSTRUCTION SEQUENCE FOR NOMINAL 
COLATITUDE MECHANIZATION 


OP Code* 

Location 

Nature of Operation 

L— - I • • •: •• • • •••••••• . Z •••••• Z I Z= HZ : ~ ~~ -J 

IAC 


Places p in Register C 

LDB 

C 3 

Loads B with (-M^L 

MUP 


x = (-M^ p) formed in A 

STA 

X 2 

Stores partial sum 

IAC 


Places x in C 

LDB 

X 2 

Places x in B 

MUP 


2 

x formed in A 

STA 

X 1 

Stores x 2 


Tdb 

X 1 

2 

Enters x in B 


MUP 

2 2 k-1 , . . A 

x - x formed m A 


IAC 


„ . 2 k-3 ... 2 k-1 . _ 

Replaces x with x in C 


LDB 

C k+2 

Enters (2 k- 1) 1 ' 1x1 B 


MUP 


Term in series formed in A 


ADD 

X 2 

Term added to partial sum 


STA 

X 2' X 

Stores partial sum or X 


— 




*The bracketed terms are repeated for each term 
in the series. 
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If y is the desired square root, then 

1 ^ y < 1. 582 

A fairly good estimate of the square root can be obtained by 


x 


o 


1 + 


s in \i 
2 


From Newton' s process, the next estimate is obtained by forming 


x 


1 


1 

2 


x + -^) 

° x o / 


Some experimentation shows that two iterations are needed to find y accurate 
to the sixth place. The sequence given in Table VI would then complete the 
calculation of \ with the great circle mechanization. 

The first command in this sequence replaces the last command of Table V. 
The memory requirements for the great circle mechanization are, there- 
fore. 


Commands: 55 + 27 = 82 
Temporary Storage: 3 
Constants: 1 0 

This sequence was formed under the assumption that the square-root com- 
mand is not part of the computer command repertoire. If the square -root 
command is available, the bracketed terms in Table VI are deleted; two 
commands replace the iteration for a net reduction of 15 commands in the 
sequence. 

The memory requirement is reduced to the following: 

Instructions: 67 
Temporary Storage: 2 
Constants: 10 

The word length of the commands depends on the number of different com- 
mands needed by the computer to conduct all of the computations. Since 
only one phase of the computer operations is considered here, one cannot 
judge command- word length on the basis of the instructions listed here. A 
fairly reasonable assumption on command-word length is 12 bits, or 
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TABLE VI. - COMMAND SEQUENCE OF, GREAT CIRCLE 
MECHANIZATION GIVEN^'SIN 


OP Code* 

Location 

Nature of Operation 

IAC 


Places sin p in C 

LDB 

C 9 

Loads Kj in B 

MUP 


Product formed in A 

STA 

X 1 

Store Kj sin p 


“SET 


Divide by 2 by Et. shift 1 bit 


ADD 

C 10 



STA 

V 

Store first estimate 

LDC 

X 1 

Load sin p in C 

LDB 

X 1 

Load sin p in B 

MUP 


2 

Form (Kj sin p) in A 

ADD 

C 10 

2 

Form Y for Newton process 


STA 

X 3 

2 

Store Y for second iteration 


LDC 

*2 

Loads first estimate in C 


DIV 


Quotient formed in B 


SLT 


Bring quotient across to A 


ADD 

X 2 

Add first estimate to quotient 


SET 


Divide by 2 


STA 


Eeplaces first estimate with second 
estimate 


IAC 


Load second estimate in C 


LDA 


Load Y 2 


DIV 


Quotient formed in B 


SLT 


Move quotient to A 


ADD 

X 2 

Add second estimate to quotient 


SRT 


Divide by 2 


_JLAC 


Y sufficiently accurate enters C 

LDA 

X 1 

Load K.^ sin p in A 

DIV 


X formed in B 

STB 

X 

Stores direction cosine 


*The bracketed terms are replaced by two other 
commands if a square root command is available. 
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one-half the length of a data word. All of the memory requirements can then 
be reduced to an equivalent number of 12-bit words; a summary of the require- 
ments are compared in Table VII. 

From this comparison, one can see that the great circle mechanization 
requires four to five times the memory required for the nominal cosine 
mechanization. If one makes a more realistic assumption (in terms of pre- 
sent-day patterning technology) that some form of polynomial mechanization 
is needed for the cosine pattern, the additional memory requirements for the 
great circle pattern are less. Current practices in microminiaturization of 
computer elements probably render the additional requirements for the 
great circle rather insignificant in terms of extra power and mass in the 
computer. 
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SECTION IV 

PICKOFF ARRANGEMENT STUDY 


An important parameter in the ESVG readout system is the location of pick- 
offs. The pickoff axis is the fundamental frame of reference of the direction 
cosine calculated from phase-angle data. For this reason, the location of 
each pickoff axis should be precisely known or precisely adjusted to its 
desired point. 

In this section, the role of pickoff arrangement to the overall contribution of 
readout error shall be studied. Two specific arrangements are selected for 
detailed study --an orthogonal arrangement of three pickoffs and a nonorthog- 
onal arrangement of four pickoffs. Since there is a desire to maintain uniform 
accuracy for all possible spin-vector directions, the pickoffs are located so 
that the distribution of axes is uniform in the three space dimensions. 

Closely related to pickoff arrangement is the readout mode; i. e. , the number 
of pickoffs employed to specify the spin-vector direction cosines with respect 
to a basic orthogonal frame of reference. When three pickoffs are employed, 
the direction cosines are completely specified; however, because the readout 
deals with a unit spin vector, the information from three pickoffs are redun- 
dant to that extent. This redundancy can then be utilized to provide partial 
compensation for any error present in the computations. When two pickoffs 
are employed, the unit vector property can be utilized to calculate the magni- 
tude of the direction cosine normal to the plane of the two pickoff axes. 
(Throughout this report this cosine is called the third direction cosine. ) The 
sign must be furnished from other information. The past history of that 
direction cosine has often been used in the past for determining its sign. 
Generally, no ambiguity develops with this technique, for when the direction 
cosine along that axis changes, another combination of two pickoffs is being 
used. 

In this section, error sensitivities for each contemplated combination of 
pickoff arrangement and readout mode are developed. These sensitivity 
equations are then combined with the pattern sensitivity equations developed 
in Section III to search for the optimum combination of parameters that 
minimizes pattern line edge uncertainty. 


SOLUTIONS TO THE THIRD DIRECTION COSINE 
SENSITIVITY PROBLEM 


The calculation of the third direction cosine from data to two orthogonally 
placed pickoffs is sensitive to errors in the results of the calculations from 
pickoff signals. It has been observed that this sensitivity varies with the 
cosines themselves, and it has produced undesirable effects in drift data col- 
lected and analyzed from present gyros. In the overall effort to improve 
readout accuracy, a study of how the arrangement of pickoffs may be used to 
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reduce the third direction cosine sensitivity and control its variation to more 
tolerable levels was conducted. 

There are several alternatives available to reduce the expected error of the 
third direction cosine. 

• Use a more precise calibration of the pattern in those lati- 
tudes where the third direction cosine error sensitivity is 
large. Such a calibration will most likely lead to a larger 
number of terms in the mechanization equation. 

• Implement a pattern wherein calculational sensitivity com- 
pensates the third direction cosine sensitivity, yielding a 
more uniform overall sensitivity. 

• Extend the pattern range so that three cosines are calcu- 
lated from pattern information for those conditions where 
third direction cosine calculational sensitivity is greater 
than desired. 

• Add and arrange pickoffs in a manner such that the third 
direction cosine sensitivity is reduced and/or is more 
uniform. 

Although the first alternative may not be the most desirable, it is the only one 
which can be employed for existing units. What experimentation that has been 
conducted has yielded improved results. It is hoped that such an alternative 
need not be further extended for future designs. The following study investi- 
gates the feasibility of the other alternatives and determines whether each or 
any combination will yield a system in which the readout sensitivity to pat- 
tern edge errors is significantly reduced. 


GENERAL DIRECTION COSINE CALCULATIONS 
FROM PICKOFF COSINES 


Let the axes of pickoffs and P 2 define a plane in the gyro stator assembly. 
Let axes and X 2 be two orthogonal axes in that plane, with axis X^^ 
located at angles and o ^ with respect to the pickoff axes as shown in 
Figure 21. 

Let the direction cosines be known, as calculated from pickoff output informa- 
tion generated by the pattern, and be denoted as cos 0 for P and cos 0 O 

for P 2 . From knowledge of cos 0^ cos 0 2 , a and Q^, it is desired 

to calculate the direction cosines of the spin vector with respect to X^^ and 

X 2 . The geometry of the calculation for the cosine with respect to X^^ is 

shown in Figure 22 . 
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With the aid of the relationships among the angles of spherical triangles, the 
equation for the cosine may be derived. From the law of cosines, 

cos = cos cos + a 2^ + sin sin (hj + d g ) cos Lj 


« j = cos cos CTj + sin 0^ sin d^^ cos 


where a ^ is the direction cosine of the spin vector with respect to X^. The 

first equation may be substituted into the second to clear it of the unknown 
angle L^. 


sino. 


a i = cos 0 


__ — r + cos 0. 

2 sin (d^ + Qg) 1 


cos d 


sin cos (d^ + o 
1 ~ sin (a^ + Og) 


Now 


sin (d.^ + a 2 ) cos - cos (a ^ + d 2 ) sin = sin (c^ + d 2 - dj) = sin a 2 


so that 




sin a ^ ‘ sin d 2 

COS 0 O T TZ T s + COS 0, — : ^ ; < 

2 sin(d 1 +CT 2 ) 1 sm (Cj+dg) 


sin d 2 


(39) 


The derivation for a^> the direction cosine with respect to X 2> is similar: 


cos d 


or 2 = cos X 2 


L 


COS d. 


sin ('d 1 + d 2 ) 


- cos 


sin (dj + d 2 ) 


(40) 


The magnitude of the third direction cosine is calculated from the orthogonality 
identity [Equation (29)]. 




( 41 ) 
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1 


( 42 ) 


a r 


= ± 


V 


1 - 


: 2 , 

Sm (CT 1 +a 2 } 


cos 0 1 + cos 0 r 

I £. 


2 cos 0^ cos #2 cos ( a j +a 2^j 


The sensitivity of a a 2 , and to errors in the calculations of cos 0^ 
and cos 0 2 is obtained by evaluating differentials of these equations. Let 


Then, 


e.^ = A cos 0.^ = - sin 0 1 A 0^ 


e 2 " A cos 0 2 - sin 0 2 A0 2 


sm a . 


A a. 


sin (Oj + (j 2 ) 


sin a 2 

e 2 + sin (a, + ct q ) £ 1 


( 43 ) 


cos a. 


A Ur 


sin (o 1 + a ) 


COS Oc 


sin (ct 1 + a 2 ) 1 


( 44 ) 


A a 


Q 9 

u 3 sin (o.^ + a 2 ) 


|cos 0 2 cos (ct 1 +-ct 2 ) - cos 0.J 


+ 


j^cos 0 1 cos (a 1 + a 2 ) " cos 0 2 


] 


( 45 ) 


ARRANGEMENT OF THREE PICKOFFS 


A uniform distribution of three axes in space is an orthogonal set of axes; the 
arrangement of three pickoffs which is considered is, therefore, an orthogonal 
arrangement. The basic reference for the direction cosines is also an orthog- 
onal set. The study of the three-pickoff arrangement is simplified if it is 
assumed that these sets of coordinates coincide. 
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Two-Pickoff Readout Mode 


The assumption of pickoff axis coincidence with the fundamental reference is 
implemented by setting 



The gyro axes direction cosines reduce to 



Qf 1 = COS Qj 

(46a) 


a 2 = cos 

(46b) 

a 3 " ± "V 

^1 - jcos 9 ^ + cos 0 2 ^j 

(46c) 

The error sensitivity equations 

are 



4o i = s i 

(47a) 


^ °2 = e 2 

(47b) 


[e x cos 8 , + £ 2 cos e 2 ] 


Thp equation for e 3 shows how its sensitivity to e 2 and ^ varies with 


spiji vector direction. If a choice q| pickoff combinations is available, the 
choice which minimize^. error ip tfee two pickoff s nearest the rotor equator, 
for then 



(48a) 



(48b) 
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If one of the two pickoffs nearest the equator is not used, one of these ratios 
is greater than the unity. 

For the comparisons that are to be made later, it is convenient to assume a 
specific two-pickoff mode. Let the spin vector be defined to be 

S = ai Q + J3j c + yk. c (49) 

A A A 

where i , j , k are the unit base vectors of the fundamental reference, 
c J c c 

Let y be the direction cosine whose magnitude is calculated from the orthog- 
onality identity. Then 


a ^ = a 

01 2 = P 

a 3 = y 


The error in the spin vector definition from readout errors is 


AS = Aa'i ( , + A/3j c + Ay k e 


./^convenient basis for comparing the error amplification is the magnitude of 
AS, which is easily calculated from its square 

AS 2 = Aor 2 + A/3 2 + Ay 2 (50) 

For the specific mode assumed, 

= A a 

5 2 = W 

e 3 = ^ 

Appropriate substitutions yield 


AC 2 A 2 

AS = A a 

\ 2 | 

1 +2—\ 
I 1 2 

+ A/3 2 

| 1 + £i| 


i y 1 


i y / 


[a <* 2 (1 - /3 2 ) + A(3 2 (1 - a 2 )] 
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The sources producing errors Aor and A/3 are such that their expected 
values are zero. Consequently, the variance of these error functions is equal 
to the expected value of thgir squares. If the variances of A a and A/3 are 
known, the variance of AS can then be calculated 

a 2 (A s) = i ^ [ a2 < X - i3 2 ) + o 2 (A|3) (1 - a 2 )] (51) 

1 - a - /3 


Three-Pickoff Readout Mode 


When the three-pickoff readout mode is used, all three direction cosines are 
determined from phase-angle information; the orthogonality identity is not 
need. In this case> the variance of AS is 

a 2 (AS) = a 2 (A a) + a 2 (A/3) + a 2 (Ay) (52) 

When this mode is used, there is no assurance that the spin vector is a unit 
vector. If it is not a unit vector, it is produced by an error in the readout. 

In this case, the orthogonality identity can be used to correct the direction 
cosines so that the spin vector is a unit vector. 

Other criteria may also be considered. For instance, the corrections may be 
weighted in proportion to the standard deviation of each of the direction co- 
sines. However, such corrections are undesirable for two reasons: 

• Weighting the corrections according to the standard deviation 
requires that the calculation of the standard deviation be pro- 
grammed in the navigation computer. One only needs to refer 
to the equations developed in Section III to see how complex 
these equations could be. Since this is likely to be the only 
place where this calculation is needed, one should be assured 
that such weighting is soundly justified. 

* Like the true spin vector, the error vector has three compo- 
nents. These components can be transformed into one com- 
ponent parallel to the spin vector and two components nor- 
mal to it. The need for q#irrection was sensed because the 
calculated spin* vector is not of unit length. This information 
has significance only on that error component parallel to the 
spin vector. 

The significance of the second reason is easily pointed out by an example. 
Suppose that corrections are weighted according to the standard deviation of 
the direction cosines in an ESVG readout system employing a great circle 
pattern on the rotor. Let the "initial" calculations from one read cycle show 
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that one of the three pickoffs is reading the pattern very close to the equator. 
From Figure 15 it is observed that the expected error from line edge disper- 
sion is greatest for that pickoff nearest the equator. However, the primary 
contribution of error from this pickoff is in the plane normal to the spin axis; 
it has little bearing on spin vector length. The spin vector is more effectively 
normalized by correcting one or both of the other two cosines, which are of 
necessity reading the extreme pattern regions in order for a three-pick<^ff 
mode to exist. Since Figure 15 indicates that the expected error from line 
edge effects for the great circle pattern are almost minimum for these two 
cosines, the implementation of weighting according to standard deviation of 
line edge errors can lead to an inconsistency between the operation and its 
purpose. 

Since the need for correction is based only on the length of the calculated spin 
vector, the most sensible correction on the spin vector is ordinary normali- 
zation: 


a 


a * 

■ \l 2 _2 

2 

(53a) 


'O' + 0 + 

y 



§ 



p' 

' v. 2 + + 

~2 

y 

(53b) 

y‘ 

V., 

“ \1 2 n 2 

~2 

(53c) 


Vq- + 0 ■ + 

y 



Such a normalization process will produce a reduction in the expected value 
_»o — 

of AS . Let the corrections based on the normalization of S be 


where 


A a* 

= a 1 - a = a | “ lj — - ■ear 

(54a) 

A/3 ' 

= 13 ' -p = fllj+s - l] * - *0 

(54 b) 

<1 

= y' -y = r (i +e - i) » - ®y 

(54c) 


a 2 + 0 2 + 2 - (1 + e) 2 

(55) 

d Ay 

are the true errors. 



(a-Aa) 2 + (0-A0) + (y-Ay) = 1 

(56) 
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If Equations (55) and (56) are combined and second-degree error terms are 
dropped, one finds that 


a Aa + /3A/3 + yA y st 6 (57) 

The error in the corrected vector is 

(AS') 2 = (Aa ' + Aa) 2 + (A/3 + A|3) 2 + (Ay' + Ay) 2 

= (Aa - ea) 2 + (A/3 - e/3) 2 + (Ay - ey) 2 

from Equation (54). Expanding and substituting (50) and (57) into the expan- 
sion yields 

(AS') 2 - AS 2 - 2e 2 + e 2 = AS 2 - e 2 
Consequently, the variance in AS ' may be found by 

a 2 (AS') = a 2 (AS) - a 2 (e) (58) 

2 

All that remains is to evaluate the variance of a (e) which can be ascertained 
from the expected value of the square of Equation (57); 

E(e 2 ) £*a 2 E(Aa 2 ) + /3 2 E(A/3 2 ) + y 2 E(Ay 2 ) + 2a/3 E(Aa A/3) 

+ 2a y E(Aa Ay) + 2/3 y E(A/3 Ay) 

The error variables Aa, A/3, Ay are correlated only to the extent that they are 
likely to be identical functions of their parent parameters a, /3, y, respectively. 
Otherwise, they are independent. Since the probability that any two of the 
direction cosines will be equal is extremely small, one can confidently assume 
that 


E (AaA/3) = E (AaAy) = E (A/3 Ay) = 0 

Consequently, 

o 2 (e) = E(e 2 ) = a 2 a 2 (Aa) + /3 2 a 2 (A/3) + y 2 a 2 (Ay) 

Therefore, 

a 2 (AS') = a 2 (Aa) (1-a 2 ) + cr 2 (A/3) (1-/3 2 ) + a 2 (Ay)(l-y 2 ) 

If there is relatively small variation of the variances of Aa, A/3, Ay as a func- 
tion of the parent paraihbters, a reasonable approximation to this equation is 

a 2 (AS') = y [a 2 (Aa) + a 2 A/3) + a 2 (Ay)] = ^ a 2 (AS) (60) 
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ARRANGEMENT OF A FOUR-PICKOFF SYSTEM 


It was noted earlier that the distribution of pickoff axes, in spatial coordinates 
to maintain accurate all-attitude readout, should be as uniform as possible. 

In the case of four pickoff s, the axes originating from the center of a regular 
tetrahedron that are normal to the faces of the tetrahedron would define the 
equal distribution of four axes in three-space. The same distribution is de- 
fined by lines which connect diametrically opposite corners of a cube. This 
fact provides the clue to the location of four pickoffs, which not only provides 
the desired distribution but also fits in well with the hexahedral electrode 
configuration - namely, at the Y-shaped intersections of the electrode boun- 
daries, A pickoff could be located at either end of a diametrically opposite 
pair of intersections. To minimize alignment stabilities, it is desirable to 
mount all pickoffs on one housing piece. The arrangement which is considered 
here is that indicated in Figure 23. 

It is worthwhile to ascertain the pattern limits required for complete attitude 
readout - one limit for which two pickoffs are the minimum number used and 
one for which three pickoffs are the minimum number used. Since the pattern 
is symmetrically located about the equator, the pattern limits will not depend 
on which end of the pickoff axis a pickoff is located. The limits can be de- 
fined with less chance of confusion by considering an equivalent arrangement 
in which a pickoff is located at each corner of one electrode. 

To assure a minimum of two pickoffs viewing the pattern, the pattern should 
extend sufficiently far that four pickoffs can view the pattern when the spin 
axis is aligned with the center of the electrode. The direction angle between 
the spin axis and each pickoff is slightly less than 55°; the pattern should then 
extend to 36° above and below the equator. The spin axis can traverse the 
path shown in Figure 24 (which lies in a plane defined by two electrode force 
axes) with two pickoffs always viewing the pattern. 

To assure that a minimum of three pickoffs view the pattern, four pickoffs 
should view the pattern when the spin axis bisects the acute angle formed by 
two pickoff axes (Figure 24). An equivalent condition occurs when pickoffs 
at diagonally opposite electrode corners view the equator; the other two pick- 
offs should also view the pattern at the same time. Either condition yields 
the same result - that the pattern should extend a minimum of 55° on either 
side of the equator. 

If pattern range has no effect on pattern accuracy and readout resolution, the 
single advantage of complete readout information (with the accompanying re- 
dundancy) from a three-pickoff readout mode would justify extending the range 
to 55° on either side of the equator. It is more likely that the pattern applica- 
tion accuracy and resolution can be increased as the pattern range is decreased. 
In this context, the treatment of pattern range as a variable in Section III is 
useful in the optimization study discussed later in this section. 
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1. Four Pickoffs View Pattern When Spin Axis 
is Here to Assure at Least Two Picicoffs. 

2. Four Pickoffs View Pattern When Spin Axis 
is Here to Assure at Least Three Pickoffs. 


Figure 24. Determination of Minimum Pattern 
Ranges for Two- and Three -Pickoff 
Modes 
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Basic Direction Cosine Equations 


The angles formed by the axes of pickoffs Pg, and Pg are nominally 

cos - * (-1/3), while the nominal angle between the axis of P' 4 with each of 

the other three is cos - * 1/3. The desired axes of reference nominally bisect 
the angles formed by the axes of P^, Pg, and Pg. Each is coplanar with, but 

orthogonal tc* the bisector of the axis of P 4 and another pickoff axis. The 

equations for force axis direction con sines should then assume a rather simple 
and symmetric character. Any computation not involving phase -angle 
data from the signals of pickoff P 4 takes the form of Equation (39) with 

o i = a g ” cos * V 1/3 

Any computation involving data from pickoff P takes the form of Equation 
(40) with 

a ^ -• Og = cos * V 2/3 


Each of the six combinations of two pickoffs from the set of four can be used 
to ascertain one of the basic cosines directly and' provides two equations for 


each direction cosine: 


_V3 
" 2 

(cos 0 1 + cos 0g) 

V 3 
" 2 

■s 

(cos + cos 0g) 

V3" 

2 

" 2 

(cos + COS 0g) 



(cos 0 4 - COS 0g) 

(61a) 

(cos 0 4 - COS 0j) 

(61b) 

(cos 0 4 - COS 0 2 ) 

(61c) 


Some relationships among the pickoff cosines are of interest and may be use- 
ful for checking the pickoff alignment. From each of these equations, one 
deduces that 


CO & © 4 — 


cosQj^ + cos 


+ cos 0g 


(62) 


The sum of the squares of the six distinct equations (two each for a, 0, and 
y) also yields an interesting result: 
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2 (u 2 + /3 2 + y 2 ) = 2 = 3/4 3 (cos 2 0^ + cos 2 9^ + cos 2 0^ + cos 2 0^) 


+ 2 (cos 9 ^ cos 9 2 + cos 0 2 cos 0^ + cos 0 ^ cos 0^) 


* 2 COS 0^ (cos 0.^ + COS 0 2 + COS 0g)j 

Using Equation (62), the following substitution can be made: 

2 cos 0^ (cos 0^ + cos 0 g + cos 0^) = (cos 0^ + cos 

+ cos 2 0^ 

so that Equation (63) can be reduced to 

O O 0 

cos 0 ^ + cos 0 2 + cos 0 3 + cos 0 4 = 4/3 

Equations (62) and (64) would be useful in aligning and calibrating the pick- 
offs. 

Two-Pickoff Readout Mode 

To illustrate the calculations needed for complete specifications of a, ft, 
and y in the two-pickoff readout mode, assume that pickoffs P^ and Pg 

are the only pickoffs viewing the pattern. Only one of Equation (61) remains 
for a direct calculation of a basic direction cosine: 

a = ~ 2 ~ (cos 0^ + 0 2 ) 

A 

The direction cosine coplanar with these pickoffs and normal to F x+ is not 
colinear with Fy + nor F z+ ; rather, it lies tt / 4 and 3tt/4 from them. 
Defining 


0 2 + cos 9^y 


(63) 


(64) 


f 12 (i y+ - F z+ ) 
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then 


F x + X F 12 


= W2 (F z+ + F y+ ) = 


34 


(65) 


The direction cosines with respect to these axes are 


S - 




(cos 


- cos 0^) 


I S • Fg^J 9/8 (cos 2 0^ + cos 2 0 2 + 2/3 cos 0^003 02) 

The equations for j3 and y are 



The sign of S • Fg^ is all that is needed for complete specification of £ 

and y . In a system where a pickoff provides information to the readout 
whenever it views the pattern, sign information can be provided by a memory 
element associated with a pickoff which is not one of the two pickoff s used 
for readout. In an orthogonally aligned three -pickoff system where the pick- 
offs are nominally coincident with the axes of reference, the application of 
such a memory element is direct. For a nonorthogonal system, its appli- 
cation is less direct, but nevertheless could be used. 


An examination of Figure 23, in view of the definition of F 34 [Equation (65)], 


shows that this axis bisects the angle formed by the axes of pickoffs 


and 


P^. In the situation where pickoffs 


and P, 


2 are the pickoffs in use, a 
memory element? which stores the sign of the last direction cosine reading 


of Pg or P^ (or both) would also define the sign of S • Fg^, That such a 
memory element can be utilized with reasonable dependability becomes 


apparent when it is observed that if the spin axis strays away from F 34 or 
-F34, either Pg or P^ begins to view the pattern, and a three -pickoff 
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readout mode or a different two-pickoff mode becomes possible. 

A A 

An alternate method of determining the sign of S • Fg 4 is to measure the 
phase angle of a pulse from generated by a meridional pattern line with 
respect to the pulse from generated by the same line. Figure 25 illus- 
trates the geometry of this method. Let that phase angle be denoted by <t> 2 y 
The sign of sin * s positive if S • Fg 4 is positive and negative if 

S • Fg^ is negative. When either one is nearly zero, Pg or P 4 or both 
will be closer to the equator than P^ or Pg, and a three -pickoff mode or 

a different two-pickoff mode is available. Thus, a measurement is needed on 
only at times when it distinctly differs from zero or rr. 

To shorten the notation, let 

„ J + 1 if sin 0n i > 0 
J 21 = \ 

- 1 if sin $ 2 i < 0 

then 
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It is readily apparent that another two-pickoff combination related to the com 
bination of and Pg is that of Pg and P^. The same direction cosine 

computed directly from Pj and Pg is computed directly from Pg and P^. 


a 



(cos 0 4 - cos 0g) 


Also, 


34 


yr 

— 4- (cos 0 4 + COS 0g ) 




6 g + cos 


a 4 - -5 cos e 3 


cos 0 4 ) 


The sign of IS’ F^g I is determined by measuring the phase of a meridional line 
pulse from P 4 with respect to the same pulse from Pg. Let 


43 


_ ) + l if sin 0 4 g > 0 
|-1 if sin 0 43 < 0 


Then 


.■4 


cos 0 4 + cos 0g + (? 4 2 


- 3(cos 2 0g + cos 2 0 4 )+ 2 cos 0g cos 0 4 


cos 0 4 + cos 0g - <? 4 g Yf - 3(cos^ 0g + cos 0 4 ) + 2 cos 0g cos 0 4 


(67a) 

(67b) 
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Equations (66) and (67) can be derived from one of the three -pickoff mode 
equations by solving Equations (62) and (64) simultaneously for one of the 
missing pickoff measurements, in terms of the two which are available and 
using the quadratic formula. A sign convention (i. e. , a 21 and a^g) still 

must be developed independently of this method for the sign of the square root 
of the discriminant. 

The sensitivity of a, /3, and y to errors in the pickoff cosines can now be 
determined. Let 


e i - A (cos 0.) 


The error sensitivity equation for a basic direction cosine can be expressed 
in general form: 



L 


2- a (cos ej 
i = 1 3 


Using this form for the combination of two pickoffs selected yields 


Vi 

A a = — gr ; (®j_ + ®2 3 


A/3 



(3 cos + cos 0 2 ) + e 2 (3 cos 0 2 + cos 



A7 



(3 cos 0^ + cos 0 2 ) + e 2 (3 cos © 2 + cos 



where 


8 2 
-g — 3(cos 0 ^ + cos 


0 2 ) - 2 cos 9^ cos 0 2 
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Under the assumptions previously made about the expectation and correlation 
of the error variables e., the variances of A a, A/3, and Ay become 


a 2 (An) = | 


^a 2 (ej) + a ' 


(fig) 


2 ,..,. 3 

a (A£S) = jg- 


o 2 (e 1 > 


1 + 


'21 


L 


3 cos&i + cos0 2 


+ a 2 (e 2 ) 


u 21 

1 ~ * — ■■ ■ 3 cos0 9 + cos0 t 

. V 


o 2 (ar) = 4 


a 2 (e 1 ) 


’211 


L ^7 


3 cos0 j + cos 0 2 


+ a> 2 ) 


a 21 

1 + _ 3 cos0 9 + cos©- 

. ^ 1 


These variances are then summed to obtain the variance of the error vector 
AS. 


a 2 (AS) = ~ 


o 2 ( £l ) 


3 + 


(3 COS0J + cos 


+ o 2 ^) 


3 + 


(3 cos0 2 + cose^ 


21 


(6 8) 
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Three -Pickoff Readout Mode 


When three piekoffs view the pattern, the basic direction cosines can be 
calculated from the three equations of (61), which do not involve the missing 
pickoff cosine. For example, if piekoffs P^, P 2 , and Pg provide signals 

to the phase -angle electronics, the central terms of (61) are used. 
Accordingly, the error sensitivity equations for this combination of piekoffs 
are 

nr- 

Ac? = ~2 (e^ + eg) 

t/sT 

Aj3 = 2~ + e 3^ 

Ay = Y (e x + e 3 ) 


The variances of the direction cosines are 


a 2 (Ac) 


a 2 (A|3 ) 


_3 

4 


3_ 

4 


a 2 (Ay) = | 


a 2 ( e ^ ) + a 2 (eg) 


a 2 (e 2 ) + o 2 (e 3 ) 


a 2 ( ei ) +a 2 (e 3 ) 


so that the variance in AS is 


a 2 (AS) = -| 


a 2 ^) + +a 2 (e 3 ) 


(69a) 

(69b) 

(69c) 

(70) 


It is interesting to note that, for the three-pickoff readout mode, AS is 
more sensitive t© a given error in a pickoff cosine in the four-pickoff arrange- 
ment than in the three-pickoff arrangement. The reason for this is easily 
noted in an expression for the sum of the squares of Equations (43) and (44): 
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1 


A 2 . A 2 

Aa 1 + A a 2 - ~—j 


sin (o 


1 
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This equation shows that the expected deviation of the direction cosines with 
respect to any two orthogonal axes in the plane defined by the pickoffs is 
minimized if the pickoff axes are orthogonal. 

The spin vector may be normalized to reduce the expected value of AS. 
Substitution for the variances of Aar, A/3, Ay from Equations (69) into 
(59) yields the desired variance for the normalized spin vector. If these 
variances are relatively constant over the range of a fairly good 

approximation is given by Equation (60). 


93 



OPTIMUM PARAMETER STUDY 


Computer Program Description 

A search for the optimum combination of pattern parameters, pickoff ar- 
rangement, and readout mode that minimizes the pattern edge dispersion 
error amplification was conducted with the aid of a computer program. The 

parameters (cos ^ a) and (cos * (3) were serially incremented over the 
pattern range being investigated; the pickoff axis direction cosines for the 
three-pickoff arrangement are calculated for each combination of values. 

The pickoff axis direction cosines for the four -pickoff arrangement are cal- 
culated with the solutions of Equations (61) for cos 0.. From the pattern 

parameters specified by the input, the expected error in the pickoff axis 
direction cosines are calculated using the equations derived in Section III. 

The expected error of the spin vector for each combination of pickoff ar- 
rangement and readout mode discussed in this section is calculated and 
printed out for easy comparison. 

The calculation of error amplification from the pattern parameters was 
coded as a subroutine to the program. In addition to specifying pattern type, 
number of pairs of lines, and the pattern range, an additional parameter was 
furnished which allows space between the extremities of adjacent pattern 
lines. Allowance for this space is needed because pattern lines must have 
a width at least as great as the diameter of the pickoff field of view and be- 
cause a finite time is required for the pulse output electronics (HIT) to be 
reinitialized in preparation for the signal produced by the next pattern line. 
The effect of the spacing margin angle was to reduce the parameter 6 in 

the error equations derived in Section III from the pattern range input to 
the program by that angle. 

Since different pattern ranges are required for all-attitude readout for each 
pickoff arrangement, provision was made to input a separate range for each. 
Ranges which were specified were equal to or greater than the minimum re- 
quired to provide a two -pickoff readout mode at all times; expected 
errors from a three-pickoff mode were calculated only if such a mode were 
possible. 


Results of Study 

It was noted on page 53 that pattern edge uncertainties produced by a 
colatitude pattern line are independent of the colatitude of the pickoff. This 
fact can be utilized to compare the amplification of error produced by the 
pickoff arrangement and readout mode. When the minimum pattern ranges 
for all-attitude readout in a two-pickoff mode are assumed for each arrange- 
ment, the combinations rank as follows in increasing order of average 
error amplification: 
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• Two-pickoff mode, four-pickoff arrangement 

• Three -pickoff mode, three -pickoff arrangement 

• Three -pickoff mode, four-pickoff arrangement 

• Two-pickoff mode, three -pickoff arrangement 
Some observations and conditions are worthy of note: 

• When three pickoffs in a four-pickoff arrangement view the 
pattern, four choices of readout modes are available : the three - 
pickoff mode and three combinations of two pickoffs. Of the 
three combinations of two, the two nearest the rotor equator 
yield the smallest error amplification. 

• When a three -pickoff mode is available, the error amplification 
in the normal spin vector from this mode, in most cases, is 
smaller than that obtained from any two-pickoff mode. This is 
generally true for the four-pickoff arrangement and invariably 
true for the three -pickoff arrangement. 

• Large variations in the error amplification occur in the two- 
pickoff modes, expecially for the three -pickoff arrangement. 

This is attributed directly to the sensitivity of the third direc- 
tion cosine error to the other cosines. The amplification is 
greatest when a three-pickoff mode is available. 

• Even when taking into account the reduced pattern range neces- 
sary for the four-pickoff arrangement, the error amplification 
from three -pickpff readout modes is smaller for the three- 
pickoff arrangement than for the four-pickoff arrangement. This 
finding is attributed directly to the fact that the pickoff axes of 
the three-pickoff arrangement are mutually orthogonal. 

Comparisons of pattern type, with all other parameters equal, showed that 
the average error amplification, in increasing order, ranked as follows: 

• Cosine Pattern 

• Colatitude Pattern 

• Great Circle Pattern 

With other parameters equal, the differences in average error between the 
cosine pattern and the colatitude pattern range from one. to six percent, 
whereas those between the great circle pattern and the colatitude pattern 
range from five to 20 percent. The reason for the significantly larger error 
amplification from the great circle pattern is not apparent. Therefore, 
the great circle pattern would be preferred only if its edge uncertainty were 
significantly smaller than that of the other pattern types. 
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Attention was focused on the comparison of the two-pickoff mode to the three- 
pickoff mode for the three-pickoff arrangement as a function of pattern 
range. It was assumed that two pairs of cosine pattern lines were used with 
a colatitude spacing margin of two degrees. The comparisons are illustrated 
in Figure 26. 

The top curve is that rms of all two-pickoff modes that were calculated; in- 
cluded were many cases for which the magnitude of either or both direction 
cosines exceeded the magnitude of the third direction cosine. The large 
variation observed in this mode is dramatized by comparing it with the 
bottom curve, which depicts the error amplification when the two pickoffs 
used for readout both view the equator. It is also dramatized by comparing 
it with the "restricted" two-pickoff rms curve, which includes only those 
cases for which | a! <;|y| and 10 | «s|y| [or when Equations (48) are 
satisfied]. Approximately 70 percent of the cases included in the top curve 
are retained in the "restricted" rms. The rms of the three -pickoff modes 
includes only the cases for which such a mode is available. As the pattern 
range is increased, the number of available three -pickoff cases becomes 
more numerous. 

The "mixed mode rms" curve is the rms of error amplification when one 
adopts a strategy of using the three -pickoff mode when it is available and 
the two-pickoff mode only when one of the pickoffs does not view the pattern. 
Note that this curve is not as sensitive to the pattern range as the others 
over the range under consideration. For the three -pickoff arrangement, 
such a strategy appears to offer the smallest overall error amplification. 

A pattern extending from a colatitude of 40° to 140° will provide a three - 
pickoff mode in all situations where a two-pickoff mode would otherwise 
yield excessive error amplification, A comparison of this strategy with 
the readout strategy employed with the present ESVG indicates a potential 
overall readout improvement of 20 percent. 

The use of two pattern lines yields minimum error amplification for all 
pattern types over the ranges which were investigated. The improvement 
it provided over one pair of pattern lines is approximately two percent for 
the colatitude line and ranges from three to five percent for the cosine line. 

If multiple pairs of pattern lines cannot be employed without significantly 
increasing the logic of the readout electronics, a single pair of pattern 
lines may well be the better choice. 

The choice of a 10° spacing margin may not be practical for all cases. 
However, since this is the angle needed to provide a usable pattern for 
the present system, it is a reasonable choice for this study. 

It was noted earlier that, with the colatitude pattern, a three-pickoff mode 
Was preferred more often than the optimum two-pickoff mode in a four- 
pickoff arrangement when the three-pickoff mode was available. A "mixed- 
mode" strategy similar to that outlined for the three-pickoff arrangement 
would then lead to the least average error amplification for this combi- 
nation. However, if the minimum range cosine pattern is used instead, the 
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optimum two-pickoff mode is found to provide less error amplification more 
often than does the three -pickoff mode. The best strategy for this combi- 
nation would therefore be one in which a two-pickoff mode is employed at all 
times; when a choice exists, the two pickoffs nearest the rotor equator must 
be determined first.. Such a strategy increases the complexity of the logic 
and leads to an average error amplification approximately two percent less 
than the mixed-mode strategy for a three -pickoff arrangement using a mini- 
mum range cosine pattern. The likely reason for such a small gain is the 
nonorthogonal arrangement employed for four pickoffs. 

It is therefore concluded that the optimum combination of parameters that 
minimizes pattern edge error amplification and retains a relatively simple 
logic in the readout electronics is: 

• An orthogonal arrangement of three pickoffs 

• One pair of cosine pattern lines 

• Pattern range extending from 4(Fto 140° of colatitude 

• Mixed-mode strategy of reading three pickoffs when 
available and reading two pickoffs only when necessary 
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SECTION V 

PATTERN EDGE IMPROVEMENT 


In the previous section we have seen how the pattern parameters and location 
of pickoffs affect the magnitude of the expected readout error from pattern 
edge uncertainties. It was shown that the realizable reduction of error when 
utilizing the pattern application technology employed in previously built gyros, 
but altering the pattern parameters, is approximately 25 percent. Such a 
reduction is insufficient to realize the potential mission requirements of the 
ESVG. Therefore, it is necessary to investigate methods of improving the 
pattern edge itself, 


SOURCES OF PATTERN EDGE UNCERTAINTIES 


A major source of edge uncertainty is provided by the pattern application pro- 
cess. The optical system detects the pattern line by the change in the intensity 
of light reflected back to the pickoff by the rotor surface. The surface on one 
side of the readout line must be highly reflective and the surface on the other 
side must diffuse and/or absorb the light. Pattern application experiments 
(which will be described in detail later in this section) have shown that greater 
accuracy in readout line definition results if the line is simple to mechanize 
and if the number of steps in the line application is kept to a minimum. A 
grit-blast method of the great circle line would, therefore, produce the least 
error. 

Another contributing source of error is the quality of the rotor surface finish. 
The polishing process tends to gouge the sintered beryllium material; conse- 
quently, satisfactory results have not been consistently obtained. Some ex- 
periments have been conducted in plating rotors with chromium before polish- 
ing. The plated rotors polished well; however, the plated surface was found 
to be unstable under large temperature variations. There is also some 
question on the stability of the plate at high rotor speeds. Nevertheless, 
further experimentation with plated rotors may provide a way to overcome 
these difficulties. 

The electric fields which support the rotor also produce stresses within the 
rotor and electrode material. If the fields become sufficiently great, clumps 
of material can break off the electrode and be deposited on the rotor. These 
deposits, commonly called arc marks, produce diffuse spots on the surface. 
These spots alter the reflectance of the surface and cause noise pulses in the 
optical signal received by the pickoff. Two aspects of arc marks are of con- 
cern in regard to readout accuracy: large marks that produce signals as 
large as that produced by a nonreflective portion of the rotor and cause false 
triggering of the signal processing electronics, and marks adjacent to the 
readout line that distort the true pulse and cause an error in the timing of 
the trigger. 

The effects of these errors, and the methods of compensation, fall into three 
categories: 


99 



• False triggering 

• Nonrandom (low frequency) edge deviation 

• Random (high frequency) edge deviation 

False triggering can.be produced by large arc marks and pock marks that are 
capable of generating relatively large noise pulses in the pickoff signal. Often 
these pulses are narrower than the signal pulses and can be filtered by pulse 
width discrimination techniques. 

Errors in the pattern line mechanization process produce deviations from the 
desired readout line. The portion that is nonrandom, or of low frequency, 
can be compensated by adding terms or adjusting the coefficients of the mech- 
anization line. It is desirable to keep these additions and adjustments to a 
minimum, as sophisticated calibration techniques would be needed to quanti- 
tively determine these additions and adjustments. 

In the grit-blast method of defining pattern lines, dispersion of the grit reduces 
the sharpness of the line edge. The slope of the pickoff output signal is reduced 
by this "fuzziness"; the triggering point as a function of colatitude becomes ran- 
domized at high frequency. A similar randomization is produced by the pres- 
ence of arc marks and polishing imperfections on the rotor surfaces. The na- 
ture of this error is such that mathematical compensation is very difficult. 
Attempts at compensating for this error have indicated that a table look-up 
method is the most practical alternative. 


RESULTS OF PREVIOUS IMPROVEMENT STUDIES 

A major area of laboratory studies for improved pattern application is that of 
vacuum or vapor depositing an optically contrasting material onto the rotor 
surface. To control precisely the region on which the material is deposited, 
a mask is needed. A deposition technique generally requires the following 
steps: 

• Apply masking material to surface 

• Remove mask from areas to be deposited 

• Clean substrate surfaces 

• Deposit material onto substrate 

• Remove mask from rotor 

The material to be deposited, in addition to being an optically contrasting 
substance, must Adhere to the surface sufficiently to withstand the forces 
from the electric fields and rapid rotation. Titanium and cerium oxide have 
been deposited on rotors and have been found to resist scuffing; whether such 
resistance indicates that the deposit can withstand centrifugal and electric 
forces is not definitely known. A method for precision- balancing an ESVG 
rotor was developed by vapor -depositing gold on the rotor. The deposit 
passed a "tape test", but was not subjected to a scuffing test; it was later 
found that this deposit was removed by the electric fields. 
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The masking material should be one that adheres to the rotor surface, is 
easy to cut away from areas to be deposited, and can easily be removed from 
the surface after deposition. For a vacuum deposition process, the mask 
should not outgas; for vapor deposition, it should not be damaged or vaporized 
by the heat. Copper has been found to be better suited for a mask material 
than other materials tested. However, it is quite fragile and can chip away 
from the pattern edge. Extreme care must be exercised at the time that the 
copper is removed from the areas to be deposited. 

The adhesive quality of the deposit has been found to depend highly on the 
cleanliness of the rotor surface and deposition system. Sufficient cleanli- 
ness was not achieved without mechanically cleaning the substrate; such 
cleaning often damaged the delicate mask. 

The overall success of this deposition technique rests on the workability and 
durability of the mask. Copper is one of the few materials which is easily 
deposited and removed from the beryllium surface; however, it appears to be 
too fragile for use as a mask. When compared with the grit- blast technique, 
the number of steps required for this technique becomes a distinct disadvan- 
tage. 

Another general area of laboratory investigation has been photo-optical film 
techniques, similar to those used in microcircuitry. However, the primary 
limitation is obtaining an adequate light source. Such a problem could per- 
haps be solved today with the use of lasers as a light source. Both the high 
output and the small divergence characteristics of gas laser beams could be 
put to good use in precisely defining a cosine line on a photo- sensitive film 
covering the rotor. Some further evaluation of these techniques would pres- 
ently be warranted. 

The presently used grit-blast technique has undergone laboratory study to 
determine how pattern line sharpness and contrast can be maximized with 
minimum disturbance of the surface. These factors depend on the size of 
the grit and the exposure to the blast. 

The size of the grit has a surprisingly great influence on the effectiveness 
of the grit- blast technique. A 10-micron grit produces the greatest contrast, 
but not a very sharp edge. A 0. 3- micron grit produces a sharp edge; the con- 
trast is not very great. A 0. 7-micron grit has virtually no effect on the sur- 
face. An investigation of the effects showed that the 10-micron grit abraded 
the surface, providing the diffuse surface that was desired. The 0. 3- micron 
grit was deposited on the surface, forming a single layer. Since such a layer 
could produce undesirable effects (such as a catastrophic arc) in an operating 
gyro, the use of this size grit was not considered further. A possible reason 
that the 0. 7-micron grit was ineffective is that this size is one that collides 
elastically with the pressed crystals that comprise the beryllium material. 

The optimum exposure time is that long enough to produce microscopic rear- 
rangement of the surface into peaks and pits. Further exposure reduces the 
contrast; this may be attributed to surface peaks being broken off by the 
excess grit. The pattern line that results becomes a depression in the sur- 
face, which can be responsible for spurious electric torques. 
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There is some speculation that nozzle geometry has some influence on line 
edge sharpness. The desired distribution of grit would indicate that a rec- 
tangular nozzle with a side parallel to the line edge would provide better edge 
sharpness than a circular noizzle. To pattern a curvilinear line such as the 
cosine line with a rectangular nozzle would require that the angle of the nozzle 
as well as its position would have to be controlled in the patterning process. 
This additional constraint was not considered to be within the scope of the 
improvement program at the time; hence, a circular nozzle is used for the 
cosine line. However, a rectangular nozzle can be utilized without great 
difficulty for a great circle pattern. 

In general, dispersion of the grit does not produce a line edge as sharp as 
could be obtained with the deposition technique. However, the grit-blast 
method does involve fewer steps and provides more consistent results than 
other techniques discussed so far. Another distinct advantage is that it 
causes the least change in rotor balance (whereas the material added in the 
deposition technique must be considered a factor of rotor balance). Conse- 
quently, it is the method adopted for patterning rotors at present. 

Other studies have been conducted to determine how the occurrence and the 
effects of arc marks can be reduced. This is an important consideration in 
the design of the suspension circuits and the electrodes; significant progress 
has been made in the long-term development of the ESVG. During normal 
operation, arcing is now likely to occur only when the gyro is subjected to 
severe acceleration environments. Little attention has been paid to the pre- 
vention of arcing when the rotor is lifted; experience has indicated that it 
deserves further study. 

In an early development program, some rotors were plated with chromium 
as a means to improve the reflectivity of the surface. The chromium did not 
gouge during the polishing process, nor did arcing damage the chrome sur- 
face as it does the beryllium surface. Chromium plated rotors offer signi- 
ficant readout improvement potential Over the unplated rotors. However, the 
plating was found to be unstable over large thermal changes because of the 
difference in the expansion coefficients between chromium and beryllium. 
Greater shape stability may be achieved with thinner coats of chromium; for 
the potential gain it offers in readout accuracy, coated rotors should be fur- 
ther studied in the future. 


ARC MARKS AND TRIGGERING ERRORS 

The pattern lihe in the presdiit system is a diffuse area on an otherwise 
reflective rotor. The height-insensitive trigger is designed to produce a 
pulse at the initiiit when the signal is a fixed proportion of the amplitude of 
the input pulse. The signal from the reflective area serves as a baseline 
for measuring the pulse amplitude. 

The presence of arc marks affects the triggering point of the height -insensi- 
tive trigger in two ways. Marks adjacent to the pattern line distort the 
critical signal on which the triggering point is defined. Marks elsewhere 
cause a change in the baseline reference, producing an error in the 
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measurement of the amplitude on which the triggering point is determined. 
If arc marks must be tolerated in the readout signal, the signal baseline 
should be less sensitive to arc marks. 

The baseline can be made insensitive to arcs if the baseline signal comes 
from a diffuse area of the rotor surface. The pulse signal must then come 
from a reflective portion of the surface. Such a system would require a 
reflective pattern line in a diffuse background. 

The diffuse background does not have to cover the entire rotor surface (this 
would not be desirable, since such a surface could encourage arcing). If 
the reflective pattern line can be rendered less vulnerable to arcs (such as 
by a chromium coat), the vulnerability of the readout accuracy to arcs and 
polishing imperfections can be reduced considerably. 


RECOMMENDATIONS FOR FUTURE STUDY 

Polishing imperfections and arc marks produce significant readout errors 
in the present system. Their occurrence or their effects must be reduced 
before accuracy improvements can be made. Some laboratory investigations 
would therefore be justified in 


• Chromium plating techniques and durability of these plates 

• Polishing techniques to improve the final rotor finish 

• Suspension modifications to avoid arcing during lift 

Some investigation of an infrared readout system would be worthwhile on the 
chance that infrared wavelengths are not diffused as much as visible wave- 
lengths by arc marks. 

Film could be deposited on the rotor and the laser could be used to define the 
line by photochemical means with a setup as shown in Figure 27. This would 
be followed by other suitable processes to complete the pattern application. 

The use of a laser as shown in Figure 27 could also be a direct means of 
pattern application if the laser has sufficient power. The beam is focussed 
so that the beams converge on the spin axis. The meridional angle inter- 
cepted by the beam remains constant as the laser traverses along z. The 
variation of beam intensity also compensates in part for the variation of the 
angle between the beam axis and the rotor surface normal. This method of 
referencing a cosine pattern has enough attractive features to warrant con- 
certed laboratory experimentation. 

Electron beams have potentially more powerful intensities. However, the 
aberrations encountered with electron lenses and space charge effects pose 
some complex problems in the control of these beams. If laser beams are 
found to be unsuitable for use in the pattern application process, electron 
beams should then be further explored as a possible alternative. 
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SECTION VI 

PHASE -ANGLE MEASUREMENT TECHNIQUES 


With the introduction of a tilted great circle line scribed on the ESVG rotor 
to provide spin-vector readout, the measurement of the meridional angle be- 
tween two points of great circle segments located at the same colatitude has 
been used to determine the location of the spin axis. To ensure adequate 
accuracy in the meridional angle measurement during the gyro test phase, an 
elaborate system has been devised for directly measuring the time interval 
between the passage of pattern lines by counting pulses of a high-frequency 
precision oscillator. Pulses generated by the pattern in the pickoff signals 
are used to start and stop the counting process. This system, in its present 
size, is not amenable for use in a vehicle navigation system. 

There is little doubt that a phase-angle measurement system based on the 
presently used scheme can be built much more compactly for use on a vehicle. 
How this can be accomplished is described briefly in this section. However, 
this high-speed, all-'digital system contains a higher degree of sophistication 
than is actually needed at the system level. Consequently, a hybrid phase- 
lock system is considered and discussed. No comparisons in physical size 
were made between the two systems in this study. The accuracy of the hy- 
brid system depends on the null sensitivity of the phase discriminator; some 
improvement is needed over the sensitivities claimed in the literature to be 
adequate for ESVG needs. Available in the hybrid system, however, is an 
analog signal that could provide rate information to the navigation system. 


ALL-DIGITAL CLOCK SYSTEM 

The phase-angle measurement system presently used in the ESVG test facility 
is one in which high-frequency precision pulses are counted; hence, the name 
"clock system". A block diagram of the system is shown in Figure 28. 

The pickoff output pulses are separated into logic pulse train P Q and phased 
trains and Pg. The phased trains are used to gate one counter which pro- 
vides a measurement of the time interval between P^ and P This time in- 
terval is proportional to the accumulated count, C^, which is then submitted 

to the computer. The logic pulse starts and stops the second counter; the 
counts accumulated on this counter, (C^ + C^)> are proportional to the rotor 

period. The ratio of to (C^ + C^) is equal to the phase angle (the dimen- 
sion of the ratio, always less than unity, is in revolutions). 

The pulse trains which gate the counters must themselves be gated by the 
readout initiation signal. The counters are first reset. Then counts are 
accumulated until M rotor revolutions are completed; the count accumula- 
tion is stopped, and the results sent to the computer. In this way, M revol- 
utions are averaged. 
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To project the physical size of a vehicle -mounted version of this system, 
three factors are considered: 

• Microminiaturization of the electronic components in terms 
of up-to-date technology 

• Instrumentation and display capability needed at the test 
phase but not in the navigation system 

• Prospect of utilizing one subsystem to provide a common 
function at several needed points in the system 

The test facility system was designed to utilize components and instruments 
available in 1962. If it were necessary to build a compact version of the same 
system for test facility use today, it is probable that a significant size re- 
duction could be effected with components and instruments now available. 

The test facility system has instrumentation and versatility not needed for a 
vehicle -mounted navigation system. Time-keeping units, display panels for 
the counters, and the circuits which adapt the unit to be used at any of a wide 
range of ESVG rotor speeds can be eliminated. Since the vehicle- mounted 
unit must be completely automated, the manual controls for readout mode and 
sampling rate would be bypassed. 


The sampling of data is initiated by a clock at precise intervals in the test 
facility system. Since the time of each sample is important in the gyro drift 
analysis phase, it is part of the data recorded with each sample. Driving 
the clock is a precision, high-frequency oscillator; the clock is essentially 
a properly scaled pulse counter and accumulator. The oscillator is needed 
to actuate the sampling for the navigation system and to provide the pulses 
needed to measure the phase angle. However, other units such as the navi- 
gation computer itself require a precision oscillator. Conceptually, it is 
possible that one precision oscillator can provide precision timing pulses, 
with proper frequency scaling for each, at all points where needed in the 
navigation system. 

When such factors are taken into consideration, a functional all-digital, 
phase-angle measurement system for a vehicle is likely at a considerable 
reduction in physical size and power from that of the test facility unit. With- 
out taking into account an oscillator for the time-keeping and time-measure- 
ment functions, the size of such a system is expected to compare with the 
size of the gyro and its suspension electronics. 


HYBRID PHASE LOCK SYSTEM 

The all-digital clock system provides phase-angle information suitable for 
immediate processing in a computer. However, the method used to provide 
this convenient form also enables a complete data sample to be taken (assum 
ing no averaging) in less time than that required for the rotor to make two 
full revolutions. The sampling rates made possible by this method far 
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exceed the system requirements, even when one averages over several con- 
secutive revolutions. Consideration of a method where comparable accuracy 
requirements may be achieved with lower- speed equipment at a sacrifice in 
the sampling rate capability has led to the suggestion of a phaselocking system 
in which the analog output of a phase discriminator is nulled with a digitally 
delayed feedback signal, A block diagram of such a system is shown in 
Figure 29. 



Figure 29. Diagram of Basic Phaselock Loop 

Pulse train P ' is delayed from P^ by time T and compared with pulse train 

Pg. The phase angle between these trains generates a voltage V' which 

through a network produces a control voltage V Q which drives a variable 

frequency oscillator (VFO). The error signal changes the frequency of the 
VFO in such a way that the delay time T is adjusted so that the phase 
discriminator produces a null signal. At null, the frequency generated by 
the VFO is a representation of the time interval between the pulse trains 
Pj and Pg. 

The delay is produced by counting pulses from the VFO. P^ starts the 

counting process. After N pulses are generated, an overflow pulse is 
produced which stops the count. The series of overflow pulses is P. 7 which 
has been delayed by 

N 

T =~~ (71) 

Under steady- state conditions, the time delay between P^ and P^ remains 
fixed. Pulse trains P^ 7 and P 2 are simultaneous, V 7 is null, and T repre- 
sents the time delay of P 2 with respect to P^ When the phase between P^ 


108 




and P2 changes, and Pg are not simultaneous, and an error signal 

(V* 1 0) is developed in the phase discriminator. The network functions 
as an integrator; a nonzero input changes V Q in the direction required to 

render P^ / simultaneous with Pg. If the loop has a relatively constant re- 
sponse time, the error signal V ' is a measure of the rate of change of the 
phase angle between P^ and Pg. 

There are two aspects of this system which require further study: 

• The accuracy to which P^ 1 can be adjusted for null output 

• The dependence of the output to a change in rotor speed 

The feasibility of using this system for readout depends on the null sensitivity 
capability to which a phase discriminator can be built. A literature search 
has been conducted to determine this capability. 


Results of the Literature Search 

Null sensitivity accuracy can be expected to depend on the input wave forms 
and the frequency range of the input signals . Literature which deals with 
phase discrimination of pulse trains at audio frequencies is considered most 
pertinent to the ESVG readout problem; however, much of the literature 
deals with sinusoidal wave forms and with sonar, radio, and microwave 
frequencies. The number of articles which provided useful information is 
small. 

A serious limitation encountered in the precise measurement of phase angles 
is the stability of the instrument supply voltages, which limits the precision 
of ordinary phase meters to 0. 1°. Kalmus and Hedrich (Ref. 2) devised a 
circuit which on alternate cycles delayed a pulse in Train A, tnen a pulse in 
Train B, by approximately 180°. A switching arrangement causes the un- 
delayed pulse to turn on a flip-flop and the delayed pulse to turn it off. The 
integrated output is a square wave superimposed on a d-c bias, with the 
amplitude of the square wave proportional to the phase angle between the 
pulse trains. In this manner, d-c shifts are eliminated from measurement 
errors. A full scale deflection sensititivity of 0.01° is claimed. The fre- 
quencies over which the authors work ranged was not mentioned. 

Of incidental interest is a circuit described by King (Ref. 3) that provides a 
polarity sensitive ternary signal for use with a phasemeter where only phase- 
angle magnitude is measured. Phase-angle magnitude is proportional to the 
pulse width of the output; the polarity of the pulse tells which pulse train 
occurs first. This information is of course vital in any phase -angle readout 
based on a null- seeking technique. 

2. Kalmus, H.P.; Hedrich, A.L. : IRE Proc. 47, 90 (1959). 

3. King, R. E. : Electronic Eng. 36, 615 (1964). 
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Accuracy claims in the literature cannot be compared directly with the ESVG 
readout requirements without taking into consideration possible errors of the 
signal- shaping circuits employed by the authors which are not within the 
scope of this measurement study. The resolution of the present ESVG pulse- 
shaping electronics is of the order of nine seconds in the meridional angle, or 
one-fourth the angle represented in the sensitivity claimed by Kalmus and 
Hedrich for their phase meter circuit. If their sensitivity includes the un- 
certainty contributed by their signal- shaping electronics, one could believe 
that their circuit may show better sensitivity when the ESVG pickoff signals 
are input. 

If no significant improvement beyond 0. 01° can be expected in the phase- 
angle measurement, one could explore methods of phase- angle multiplication 
which would render an effective phase meter sensitivity commensurate with 
the resolution of the input signal. A possible method is suggested by the 
block diagram of Figure 30. 
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Figure 30. Phase-Angle Amplification by Doubling the Frequency 


The input signals are delayed 180° and superimposed onto the original signal, 
which doubles the pulse frequency. On the output, the pulse period is 

one-half that of the input period P^_ while the delay time T is maintained. 

The phase angle, proportional to T/P^, is doubled. Delays of this nature 

could be applied successively as long as T > P, , A "coarse-fine" measure- 

K I * 

ment can be obtained by using the phase meter output from A K and B K as 
the fine measurement and the pulse count in A k before a pulse occurs in B° 

k o 

(or pulse count from B before a pulse in A ) as the coarse measurement. 

The phase nulling process by the feedback circuit would be a two-step process, 
a nulling of the coarse measurement followed by a nulling of the phase meter 
output, 
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Compensation for Rotor Speed Changes 
The phase angle is related to the delay time at null as follows: 

0 = 2tt f T 

where f^ is the rotor frequency. Substituting this relation into (71) yields 

2tt f N 

V = — f L - (72) 

To mechanize a measurement for v, a provision to evaluate the ratio f^/f 

is needed. A pulse counter similar to the units used in the clock system could 
be used. Pulses which are a multiple of rotor frequency may be accumulated 
in a counter which is gated and reset by pulses from the variable frequency 
oscillator. For the required resolution, the multiple of f to be counted 

5 r 

should be of the order of 5 x 10 times the range of f . An alternative would 
be to divide down f with a countdown similar to that in the feedback loop in 
Figure 29. 


An alternative is to allow N to be a rotor speed dependent variable. 
If 



then 


v 


2ttk 

f 


(73) 


and f is a measurement of the reciprocal of v. The reciprical result, in 
the case of a cosine pattern, is not the most desirable, but this can be recti- 
fied in the mechanization equation for the direction cosine. For a great circle 
pattern, the cosine mechanization from the reciprocal of phase is not likely 
to be any more complicated than that from the phase angle itself. An unat- 
tractive feature of this alternative is the calculation of N from f with (73); 
this would require a division process in the computer, or countin§ pulses 
from a reference oscillator with gating provided by reference pulses from 
the rotor. 


The most practical method of calculating v with (72) is to mechanize that 
equation in the computer. Low- speed counters can be used to measure f 
and f, or parameters that enable their quick calculation. 
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A Phaselock System for ESVG Readout 

A block diagram of a potentially feasible phaselock system with compensation 
for rotor speed is shown in Figure 31, Signals containing logic pulse P Q 
and phased pulses P^ and P 2 are separated; the logic pulse is used as 

the basis of separation and for providing rotor speed information. The same 
logic pulse can be used with the logic pulse from another pickoff to determine 
the sign of the third direction cosine, in the manner discussed on page 87 
with the King polarity sensitive circuit. 

To provide adequate feedback resolution, the VFO output frequencies are 
likely to be quite high. It may be found convenient to heterodyne this output 
f down to a lower frequency f '. This output and the logic pulse frequency 
multiplier output would trigger pulses which are counted and input to the 
computer for calculation of the phase angle by Equation (72). The precision 
clock pulses can be input under computer control and be used to start the 
readout cycle. 

In addition to driving a network that adjusts the VFO output, V' is fed into a 
second network which outputs a signal proportional to the phase- angle rate 
of change. The nominal purpose of this network is to compensate for any 
feedback loop nonlinearities. 

The successful implementation of a phaselock system depends primarily on 
the null sensitivity of the phase discriminator. With the implementation of 
the circuit reported by Kalmus and Hedrich, there is a reasonable chance 
that a system compatible with accuracy requirements and/or other sources 
of error can be built. If it is insufficiently accurate, a phase amplification 
circuit (by stepping up the pulse frequency as shown in Figure 30) can be 
considered. To be sufficiently accurate, high-frequency pulses must be 
counted in the countdown delays; frequencies of the same order as the VFO 
output would be needed. The phase-angle nulling process would consist of 
a two-step "coarse-fine” sequence. The phaselock-loop block diagram would 
be similar to that shown in Figure 32. 


CONCLUSIONS 


Before an experimental phaselock system can be built and tested, further 
study is needed to determine the loop and network parameters. If the simpler 
loop of Figure 31 can be implemented with sufficient accuracy, chances are 
that the physical size of a phaselock system can be significantly smaller than 
the all-digital clock system. If it is important to obtain rate information 
from the ESVG, and a phaselock system indeed provides this information, 
such a system is certainly an attractive one. 

The additional hardware shown in Figure 32 would increase the physical size, 
how much depends on the size requirements of the frequency multipliers. 
Under these circumstances, other criteria than physical size may determine 
the better choice between the clock and the phaselock systems. 
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COMPUTER 



Figure 31. A Phaselock System for ESVG 
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Figure 32 . A Phaselock System with Phase Amplification 












SECTION VII 

DYNAMIC ENVIRONMENTS 


Maximum readout precision is possible only if the ESVG rotor maintains the 
position in the electrode cavity that it assumes during the readout alignment 
calibrations and adjustments, and if the pickoff colatitude does not change 
during the readout cycle. In this section, the effects of environments which 
produce rotor displacements and changing colatitudes are evaluated - static 
acceleration, vibration, and angular rotation. These effects are minimized 
if readily mechanizable computer compensation for translation is implemented 
if several rotor revolutions are averaged per data cycle, and if all 
pickoff s are read out simultaneously. 


STATIC ACCELERATION 


The ESVG rotor translates with respect to the cavity center in response to 
vehicular acceleration. This translation will cause a change in the colatitude 
which some pickoffs of the readout system view. 


Geometrical Considerations 

A rotor whose center, O', is displaced from its nominal location O in 
response to an acceleration vector X is diagrammed in Figure 33. With 
respect to a spin axis-oriented cartesian coordinate system center at O', 
the direction of pickoff P 1 has changed. On the other hand, the direction 

of pickoff Pg» which is in the direction of A, has not changed. Where P^ 

may be viewing a different rotor colatitutde, the only change occurring to Pg 
is in the distance between the pickoff and the rotor surface. While this has a 
caused a change of focus, the change is small enough (a few microinches 
compared to a distance of 0. 1 inch from the pickoff lens to the rotor surface) 
that its effect is not noticed. Therefore, the component of displacement in 
the direction of the pickoff axis causes no error in the colatitude reading. 
Only the projection of displacement in the plane normal to the pickoff axis 
needs to be considered. An equation for this projection is 


— * 
T. 

l 


A A 
(P. x T) x P. 

l l 


(74) 


where is the translation vector and P^ is the unit pickoff axis vector. 

Since it depends on the pickoff axis in question, the projection is indexed 
accordingly. 
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Figure 33. Geometry of Rotor Response to Acceleration 

The change in colatitude observed by the pickoff depends also on the orienta- 
tion of the projected translation with respect to the spin vector. For example, 
if the spin vector is normal to the plane of Figure 33, no change in colatitude 
would be sensed by pickoff P^„ It is conventient to consider components of 

T. with respect to coordinates defined from the spin vector; for this purpose 

the projection of the spin vector in the plane normal to P i is also defined: 

3. = (P^xS) x P. (75) 


As shown in Figure 34, T. can be resolved into components parallel to and 
-» i 

normal to S.. Before considering each of these components separately, it is 
1 A A 

worth mentioning that when S is parallel to P^, S. is zero. This case need 
not be separately considered because pickoff P^ is not being used, since it is 


not viewing the pattern. Consequently, will never be zero for realizable 
readout conditions. L 


Readout Translation- Sensitive Error Equations 

-4 *4 

The parallel component of T^ to (denoted by p ) changes the observed 

colatitude by the magnitude of p, as measured on the surface of the rotor, 
and is independent of the magnitude of §. . Consequently, this component can 

be evaluated by computing the inner product of T. and S.» where S. is the 

unit vector parallel to From Equations (74) and (75), 
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Figure 34, Observed Colatitude Sensitivity to Rotor Translation 


p .= T T~ [(P, x T) x P.] • [<P. x S) x P.] = r A0; 

| P.xS | 1 ii i 

Let 0 ^ denote the eolatitude read by P^ if the rotor were centered, or at 
zero displacement. Then 

A A 

|P. x S I = sin 0. 

' i ' i 

Using some vector product identities, 

[(P. x f) x P.] • {(P. x S) x P.] = [P. x (T x P.)].* [P. x (S x P.)] 

= C(T • P.) P. - T ] • [(S • P.) P. - S] = T • S - (T • P t ) (S * P t ) 

= T • S - cos 6. (T • P.) 

l i 

A A 

where S * P^ = cos 0^ , 

If the direction cosines are the basic spin vector variables rather than the 
colatitude angles, the change in direction cosine is a more desirable entity: 

1 _ , A A 

A(cos 0.) = sin 0. A0 t = — T * (P. cos 0. - S) (76) 

Note that this change is expressed in terms of known or measurable param- 
eters. The translation vector T can be determined from capacitance bridge 
measurements in the suspension, or it can be evaluated from the equation 
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T = G A 


where G is the suspension gain matrix and A is the acceleration of the 
gyro. For an isotropic support system, G reduces to a scalar. The accel- 
eration vector can be measured with accelerometers. As a result, mathe- 
matical compensation for direction cosine error from rotor translation pro- 
duced by steady-state acceleration can be implemented. 


A potential difficulty with using Equation (76) for compensation is that it is 
written in terms of the spin vector coordinates with the rotor centered, in- 
stead of the observed coordinates. However, if the change is small enough, 
the observed coordinates can be substituted for the centered spin vector 
coordinates without incurring intolerable error. From Figure j) 3 it is seen 
that the magnitude of the change does not exceed the ratio of |T | to r. In 
Volume HI, the maximum displacement of the rotor is assumed to be 10 
percent of the rotor-electrode gap; in Volume II it was observed that the 
ratio of gap to radius on existing gyros ranged from 1:167 to_I;500. If one 
makes the rather pessimestic assumption that the ratio of |T j to r does 
not exceed 10"^, the error of the compensation would be equal to the square 

“6 

of that ratio, or not greater than 10 . This is equivalent to an uncertainty 

of a fraction of an arc second and can, therefore, be neglected. 


The normal component of T\ to (denoted by n) produces a change in 

pickoff colatitude resulting from the curvilinear nature of colatitude paral- 
lels on the rotor surface. Referring to Figure 35, the change in direction 
cosine can be expressed as 



Figure 35. Curvilinear Shift Geometry 
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Therefore, 


2 

A (cos 0.) = COS 0.(COS — - 1) =“ - n -7T COS 0. 

1 1 r 2r 2 1 

Similar to the parallel component p, the normal component n never exceeds 

— * * 

| T | . Under the same assumptions about the relative magnitude of j T | to r, 
one concludes that compensation for the curvilinear effect, while it is possible, 
is not needed since the error is small enough to be tolerated. 


VIBRATION 

In a vibration environment, the acceleration is a bounded function of time. 

The translation function determined by the suspension servo transfer char- 
acteristics is also a bounded function of time. The equation for tjje direction 
cosine change from translation (76) can then be considered with_T as an in- 
dependent variable in time. In inertial space, the spin vector S remains 
invariant except for its drift from torque. Over the interval of |ime required 
to read the gyro, the drift can be neglected. The direction of will change 
only if rotation or rotational vibration is present.^ Hence, if translational 
vibration is the only form of vibration present, P. and the true colatitude 
angle 0. are invariant, and the apparent change indirection cosine is the 

inner product of a constant vector and T. Since the direction cosine variation 
is a scalar, it can be boiled down to a Scalar function of the vibration: 

A(cos 0^) = f(t) 

Since no products of vibration-induced variables occur, the characteristic 
frequencies of the direction cosine variation are the same as the vibration - 
induced translation. 

The observed colatitude angle, however, will change with vibration; for that 
reason, it is worthwhile to evaluate the effect of using observables in the 
compensation scheme for vibration inputs. Such effects depend on the vibra- 
tion frequency. Three different ranges are considered: 

• Frequencies that are low compared to the frequency repre- 
sented by the data cycle 

• Frequencies above the response frequency of the translation 
measurement apparatus 

• Intermediate frequencies 

In this discussion, it is assumed that the response frequency of the transla- 
tion measurement apparatus is higher than that represented by the data cycle 
(equal to the reciprocal of the data period). 
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At the low frequencies, the translation and colatitude change is small during 
the data period. The steady- state component of both parameters is expected 
to be of greater significance in the correction. As such, the compensation 
error is of the order of that of steady-state acceleration, which was previously 
shown to be less than 10"®. 


At the high frequencies, the translation measurement apparatus doesn't follow 
the vibration; high frequencies are then not expected to be sensed by the out- 
put. Consequently, no correction for high-frequency vibration is possible, 
and one must depend on the spin vector observables to assume an average 
value of its instantaneous function. Except at frequencies resonant with rotor 
frequency and its harmonics, some averaging can be expected. The standard 
deviation from the mean of a sample of data points taken during high-frequency 
vibration will likely be greater than that of a similar sample taken during no 
vibration. 

At intermediate frequencies, the compensation effects are not as predictable 
as for the other regions. Errors larger than steady-state errors can occur; 
for instance, if the spin vector observables are some average values of the 
instantaneous values and if the translation measurement yields an (instantan- 
eous) extreme value during a data sample. If such is indeed a possibility, it 
may be desirable to reduce the response frequency of the translation measure- 
ment apparatus to one approximately equal to the data cycle or to bypass the 
computer compensation for vibrations above a specified frequency. 

A 

When rotational vibration occurs, and cos 0^ are functions of the vibra- 

tion, and their product contains rectified terms that are second degree in the 
vibration amplitude. Since rectification occurs, the average value of the 
product is not zero, but deviates from zero in proportion to the square of the 
vibration amplitude. If the amplitude is small, such products can be disre- 
garded. 

If rotational vibration occurs without translation vibration or steady -state 
vibration, the translation vector T is zero and no correction is available. 

An increase in the standard deviation of a sample taken during this vibration 
can be expected. If either translational vibration or steady acceleration 
occurs with rotational vibration, a rectification occurs which is at least 
second degree in the acceleration and vibration parameters. The specific 
vibration parameters must be known before the rectification can be realist- 
ically assessed. Generally, these effects are small enough to be tolerated 
without compensation. 
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ANGULAR ROTATION 


When the gyro case is rotating, a pickoff direction cosine is constantly chang- 
ing. When the direction cosines are read out, resolution is lost, with the cal- 
culated cosine being some value which occurs during the readout interval. An 
upper bound of the change of colatitude observed by a pickoff during a readout 
interval is 


where Q is the vehicular rate and t is readout interval. Generally, it is 
small enough that over the change in colatitude the change in phase angle is 
linear with time: the measured colatitude from the pickoff data, therefore, is 
expected to be the colatitude at the midpoint of the data cycle. If the data 
cycle mid points for all pickoffs can be made to occur simultaneously, maxi- 
mum precision in specifying the spin vector direction during rotation is 
assured. 


Equations of Motion 

Consider a rotating set of coordinates in inertial space. By the coriolis law, 


where 


, a d / A 

Tt s m -3T s+QxS 


_d_ 

dt 

d ' 
dt 


Q 


time derivative with respect to the inertial frame 

time derivative with respect to the rotating frame 

angular velocity of the case frame with respect to the inertial 
frame 


It is assumed that the unit spin vector is space stabilized. 

Thus, 

■ »• 

and the observed rotation of the spin vector in the rotating frame is described 

by 

-Jf-s = S x n (77) 

dt 


Let 
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s = a i c + 0 a c + y k c 

AAA 

^ = n i l c + h + n 3 k c 

where a, 0, y are the direction cosines of the spin vector in the (rotating) 
case coordinates. 


Then 


4f<“ ‘ G + 13 j c + y k c> ■■ <“ 4 c + P j c + y k c> x 4 c + ‘Vc * n 3 k c J 


Now 


because 


A! g _ A!°L \ + silfL ] + AlX. Z 

dt dt c dt dt K c 


4i “i -44-j * 44 i = 0 

dt c dt J c dt c 


which is another way of saying that the rotating coordinate unit vectors are 
not moving in the rotating system. 

The vector equation reduced to three scalars: 


^ ,Qi st So 
dt pu 3 

- yfig ~ & 

(78a) 

^ = yn i 

- aQ 3 = 0 

(78b) 


- 0^ = y 

(78c) 


Now, information on the spin vector components can be provided by the read- 
out system and the applicable mechanization equation. If the phase-angle 
measurement electronics utilize a phaselocking technique like that described 
in Section VI, information on the rate of change of a , 0, and 7 may be ob- 
tained rather simply. The network driving the variable frequency oscillator 
is basically an integrator of V 7 with respect to time; consequently, V 7 is a 
function of the rate of change of the phase angle. If the relationship between 
V' and phase-angle rate is linear, and a cosine pattern is used, V 7 only 
needs to be properly scaled to provide a direct output for a , 0, or y . 
Consequently, information can readily be made available for the solution of 
Equations (78) for the angular rotation. 
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Gyro Information Limitations 

Complete specification of Q from the solution of these equations is, however, 
impossible. Although a, (3, and y may be independently determined with a 
three-pickoff mode, they are not independent because they satisfy the identity 

a 2 + p 2 + y 2 = 1 

The derivative of this identity with respect to time is 

aa + (3(3 + yy = S • — — S a 0 

This is merely an equational statement of A the fact that no information regard- 
ing changes in rotation angle parallel to S is available from the gyro; a 
second gyro whose spin vector is nominally orthogonal to § is needed to 
provide that information. 

— > 

That complete information about Q is not possible from a solution to Equa- 
tions (78 ) is readily apparent because Equation C/7)^ from which the scalar 

system is derived, deals with the cross product SxQ. Information about 

S * Q is immediately lost in this operation. If complete information about 

Sx Q is obtained from a solution of (78), the projection of Q in the plane 

A 

normal to S can be obtained from the triple product 

(S x q ) x § 

When similar information from a second gyro (whose^spin vector is nominally 
orthogonal to S) is added, a complete solution for n is, in principle, 
obtainable. 

A solution of the system (78) will not be attempted here; some characteristics 
of a, (3, and y with the time will be derived if Q is constant. 


Constant Angular Rotation 

— * 

Assuming that Q is constant. Equations (78b) and (78c) can be substituted 
into the derivative of (78a) with respect to time to give 

a + a(Q 3 2 + Qg 2 ) = (y Qg + jBiy (79) 


Now, if 

Q • S = 0 (80) 


then 

yQ 3 + 0C^ = -a: ^ 
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0 


( 81 ) 


and a differential equation in a only is obtained: 

a + a( Q 3 2 + n 2 2 + Qj 2 ) = 

The same manipulations give similar differential equations for £ and y. 
Thus, the direction cosines are sample harmonic in time if the angular ro- 
tation is constant and normal to S. Equation (79) can now be differentiated 
with respect to time and (78b), (78c), and (78a) in sequence substituted into 
that derivative to obtain the same differential equation in or; 

a + a( Cl^ + Cl 2 ^ + = 0 (82) 

Consequently, the first derivative of the direction cosines satisfies the equa- 
tion for simple harmonic motion. The general solution to a is, therefore, 

a = a + A q cos (Cl t - xfy) 

where 

o = i Hj - 


It is interesting to note that 0 can be determined by measuring the period of 
the direction cosine function. From this calculation and a solution for 

A A — ► 

(S x Cl), the magnitude, but not the sign, of S • Cl can be determined. If the 
apparent change of rotor speed due to the angular rotation can be sensed, the 

sign of the change would determine the sign of S ' Cl . Neither of these deter- 

minations would be precise, except in the case that |S : Cl | is large com- 

A 

pared to |S x Cl | . 

In examining the general solution (82)_in yiew of (79), (80), and (81), it is 
apparent that a is non zero only if Q • S is non zero. From the geometry 
of conditions imposed by constant fi, one intuitively arrives at the following 
equation for a: 

a =- Q ir (§ • Q) ( i c • Q) 
similar equations for ^ amd y can also be written. 


CONCLUSIONS 

From this study of dynamic environments, some general conclusions can be 
stated; 

• Static acceleration produces a rotor displacement which 
produces an error in readout. However, by utilizing the 
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parameters ordinarily measured in a navigation system, 
adequate mathematical compensation can be implemented 
with a computer. 

• Vibration causes a rotor displacement that varies with 
time within bounded limits. This causes a variation in 
phase-angle measurement that can be reduced by aver- 
aging over several revolutions and/or over several 
samples. Even so, the dispersion .of readout samples 
from their means is likely to increase with vibration. 

• Angular rotation produces a constantly changing direction 
cosine. When the rotation is constant, the direction cos- 
ines are sinusoidal function of time. To produce reason- 
ably precise readout of the spin axis, the pickoffs must 
be read out simultaneously. If direction cosine rates are 
measured, information ordinarily provided by a rate gyro 
can be calculated from ESVG readout parameters. 


125 



